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Chapter 1 
Introduction 



The Standard Model (SM) P is very successful in describing the known phenomena 
of particle physics. However, it also has some drawbacks. Supersymmetry (SUSY) [2[|3], 
which is a symmetry which connects fermions with bosons, is one of the best motivated 
extensions of the SM. 

In the Minimal Supersymmetric Standard Model (MSSM) [2|,[3], we introduce a super- 
partner to every known particle, i.e. scalar superpartners (sfermions) to the SM fermions, 
fermionic superpartners (gauginos and higgsinos) to the gauge bosons and Higgs bosons. 
After electroweak symmetry is broken, fields with same quantum numbers can mix. We 
therefore have sfermion mixing, and higgsino and gaugino mixing into charginos and 
neutralinos. Moreover two Higgs doublets are necessary in the MSSM. This leads to five 
physical Higgs bosons after the electroweak symmetry is broken. If i?-parity is conserved, 
the lightest supersymmetric particle (LSP), which in many scenarios is the lightest neu- 
tralino x?) appears at the end of the decay chain of each supersymmetric particle. The 
LSP escapes the detector, giving the characteristic SUSY signature of missing energy. If 
SUSY exists at the electroweak scale, experiments at future high energy colliders should 
be able to discover the superpartners of the known particles, and to study their properties 

At the Large Hadron Collider (LHC), the total SUSY production-cross section is 
expected to be dominated by the production of gluinos and squarks, which decay into 
lighter charginos or neutralinos. Of particular interest are decay chains leading to the 
next-to-lightest neutralino X2- X2 turn can always decay into the LSP Xi two 
fermions //, at least for light SM fermions. Depending on the neutralino, sfermion and 
Higgs boson masses, the possible decays of X2 three-body decays X2 ~^ Xif f 1 cascade 
two-body decays xl^ ff ^ Xiff and/or X2 ^ ^ X?//, where 0° stands for 

one of the neutral Higgs bosons or the neutral Goldstone boson of the MSSM. The leptonic 
final states are of particular interest, since they can be identified relatively easily even at 
the LHC. Moreover, the dilepton invariant mass distribution can be measured accurately 
at the LHC. In particular, its endpoint is used in several analyses which aim to reconstruct 
the mass differences of the supersymmetric particles [HE]. Under favorable circumstances 
it has been shown that this endpoint can be measured to an accuracy of 0.1% at the 
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LHC [4j. In order to match this accuracy in the theoretical prediction, at least one- loop 
corrections to X2 decays have to be included. At the planned e^e~ linear collider ILC |5], 
the lighter supersymmetric particles can be produced directly. The detailed analysis of 
X2 decays can then yield information about the supersymmetric particles. Moreover, the 
masses of the supersymmetric particles are expected to be determined with high precision 
at the ILC [5], again making the inclusion of quantum corrections mandatory to match the 
experimental precision. 

Leptonic two-body decays X2 ~^ ^f^^ Xi^^^^ have been investigated at tree level 
in Ref. [B], where /i stands for the lighter slepton of a given flavor. Three-body decays 
of X2 have also been studied at tree level in Refs.p |10]. In this thesis, we calculate 
leptonic X2 decays at one-loop level. Cases where two-body decays X2 ~^ ^1^^ ~^ Xi^~l~^ 
are kinematically alowed are treated both completely and in a single-pole approximation. 
In the complete calculation one has to employ complex slepton masses in the relevant 
propagators and one-loop integrals. The single-pole approximation in this case is performed 
in the way that the X2 decays are treated as the production and decay of the sleptons 
li. We compare the results from the complete and approximate calculations and find a 
good agreement. We also analyze a scenario where X2 only has three-body decays. The 
calculations for these decays are similar to the complete calculation except that we do not 
have to introduce complex masses. The differential decay width of X2 ^ function of the 
dilepton invariant mass Mi+i^ is calculated. The shape of the Mi+i- distribution is expected 
to be altered by real photon emission contributions, which must be added to the one-loop 
corrrections in order to cancel the infrared divergences in the virtual contributions. In 
order to obtain the total decay width of X2 ^^d hence the branching ratios of its leptonic 
decays, the invisible decays X2 ~^ Xi^i^i ^^^^ the hadronic decays X2 ~^ XiQQ ^-^e also 
calculated. 

The general MSSM has more than one hundred unknown free parameters. We assume 
CP-conserving MSSM with real parameters. For specific numerical evaluation we consider 
the SPSla parameter set [Zl IB], which is one of the standardized benchmark scenarios. 
It gives rise to a particle spectrum where many states are accessible both at the LHC 
and at ILC [6j. Note in particular that the two-body decays X2 ~^ ^f^^ ~^ Xi^^^^ 
kinematically allowed in the SPSla scenario. No other two-body decay mode is open. 
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The outline of this thesis is as follows: 

In Chapter [2] the basic ideas of SUSY are presented first, where we introduce the SUSY 
algebra, the superfields and the SUSY Lagrangian in detail. Then the MSSM Lagrangian 
and the soft SUSY-breaking Lagrangian are constructed. At last all the physical fields of 
the MSSM are discussed and the expressions for their masses are presented. 

The MSSM is renormalized in Chapter [3] where we follow the strategy of Refs. fHX . 
All relevant parameters in the MSSM are assumed to be real quantities. This amounts 
to the assumption that the soft SUSY-breaking terms conserve CP. We introduce two 
renormalization schemes, the DR and the on-shell scheme. The SM sector, the chargino and 
neutralino sector and the sfermion sector are renormalized via the on-shell renormalization 
scheme, while the Higgs sector is renormalized via a mixing of on-shell and DR schemes. 

In Chapter mthe leptonic decays X2 ^ Xi^^^^ are calculated at one- loop level. The tree- 
level calculations for these leptonic decays are outlined in Section 14. ![ where we calculate 
the leptonic decays in both a complete and an approximate way where li can be on shell. 
In Section S2] we discuss how to calculate these decays completely at one- loop level, where 
the virtual corrections and the real photon bremsstrahlung are investigated in detail. The 
approximate one-loop calculations for these decays where /i can be on shell are given in 
Section 14. 3[ We study the total decay width of X2 the branching ratios of the leptonic 
decays in Section 1^7^ In Section the explicit values of the MSSM parameters and the 
spectrum of the supersymmetric particles at the SPSla benchmark point are given first. 
Then the numerical results of the calculations are presented. 
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The Minimal Supersymmetric 
Standard Model (MSSM) 

2.1 The Standard Model 

The Standard Model (SM) of particle physics describes the electroweak and strong interac- 
tions with very good precision. It is based on the gauge group SU{3)c x SU (2)2, x f/(l)y [1] , 
where SU{3)c and SU{2)l x U{1)y describe the strong and electroweak interactions, re- 
spectively. 

The elementary particles include spin-^ fermions and spin-1 gauge bosons. The interac- 
tions in the nature are described as the exchange of different gauge bosons. The SM gauge 
bosons are eight gluons 6*^(8, 1, 0), three weak bosons W^{1, 3, 0) and a hypercharge boson 
-8^(1,1,0), where we have indicated their quantum numbers with respect to the gauge 
group SU{?>)c X SU{2)l x f/(l)y. These gauge bosons are the mediators of the strong, 
weak and electromagnetic interactions. 

The SM fermions are chiral because the left-handed and right-handed fermions trans- 
form as doublets and singlets of the gauge group SU{2)i, respectively. The fermions, i.e. 
leptons and quarks, can be arranged in three generations as below, 

1st generation 2nd generation 3th generation group representation 
eR, fiR, tr, (LL-2) 



u \ I c \ I t 



1 



(3,2,1/3) 



d I \ s J / \ b , , 

/Li \ / Li \ / Lj 

ur, cr, tR, (3M,4/3) 

dR, sr, bR, (3M, -2/3) 



In the SM we introduce one scalar doublet H with hypercharge Y = +1, which imple- 
ments the spontaneous breaking of the electroweak SU{2)l x f/(l)y symmetry down to 
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the electromagnetic U{1)em symmetry, 



SU{3)c X SU{2)l X f/(l) 



SUi3)cy<U{l)EM. 



This is called the Higgs mechanism ^14j. which gives masses to the fermions and produces 
three massive vector bosons and Z, and a massless photon 7. The Higgs mechanism 
also predicts a new particle: the Higgs boson. 

The SM provides an extremely successful description of the known phenomena of par- 
ticle physics. All the particles that it predicts, except for the Higgs boson, have been 
discovered experimentally. The mass of the Higgs boson is severely constrained from elec- 
troweak precision data[15]. Moreover, the present experimental data agree with the SM 
predictions very well|15[ [16] . However, the SM also has some drawbacks and unsolved 
problems. Here we focus on why we need supersymmetry. 

The SM is an "effective low-energy theory" for energy scales up to 100 GeV. More 
reasonably, we regard the SM as part of a larger fundamental theory which describes 
physics at arbitrarily high energies. The Quantum gravitational effects become important 
at the Planck scale Mp ~ 2.4 x lO^^GeV. A Grand Unification Theory (GUT) may appear 
at a somewhat lower energy scale. The surprising thing is that the ratio ~ 10^^ is so 
huge. This is called "hierarchy problem" |17j. 

Moreover, radiative corrections must be included in the theoretical predictions in or- 
der to match the precision of experimental measurements. When we consider radiative 
corrections to the Higgs boson mass, mfj is found to be quadratically divergent at one- 
loop level[T8]. For example, consider the contributions to m'jj from a loop containing the 
fermion / as shown in Figure 2.1(a) and suppose the Lagrangian term is —XfHff. Then 
the correction is given by 



where A is a momentum cutoff used to regulate the loop integral, which represents the 
energy scale at which new physics appears. If A is of order of the Planck scale Mp, the one- 
loop quantum correction Am^ is some 30 orders of magnitude larger than 171% < (ITeV)"^, 



Ami 



H 




A^ + logarithmically divergent terms , 



(2.1) 




s 





Figure 2.1: Radiative corrections to the Higgs boson mass 
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which is needed to preserve unitary in the scattering of the longitudinal gauge bosons [19]. 
This is the technical aspect of the hierarchy problem. One needs extreme fine-tuning of 
the parameters to cancel the large quadratic contributions against a counterterm, leaving a 
resulting Higgs boson mass of about 1 TeV. Such a cancellation is unnatural and in general 
not acceptable. Additionally, in every order of the perturbative theory, the parameters 
must be re-tuned. 

The only known way to cancel these quadratic divergences is by introducing a new 
partner with spin differing by | for every known particle. The properties of the known 
particles and their new partners are related by a symmetry, which is known as supersym- 
metry. As shown in Figure 2.1(b), suppose there is a new partner 5* for fermion /, which 
couples to the Higgs boson with a Lagrangian term —XsH'^S'^, then the contributions are 

Arrijr ~ — ^A^ + logarithmically divergent terms. (2.2) 

If each of the SM quarks and leptons has a new supersymmetric partner with A5 = X'j, 
then the quadratic divergences neatly cancel. 



2.2 Super symmetry 

In order to protect the Higgs boson mass from quadratically divergent radiative corrections, 
we introduce supersymmetry(SUSY), which is a symmetry which connects fermions and 
bosons. In this section we first introduce the SUSY algebra and the different superfields, 
and then come to the construction of the SUSY Lagrangian. The breaking of SUSY is also 
discussed. The discussions in this section are based on Refs. [21 [3l [20] . 



2.2.1 SUSY Algebra 

The generators of SUSY must turn a fermionic state into a bosonic state, and vice versa, 

Q\boson) = Ifermion), Q\ fermion) = \boson). 

This implies that the SUSY generators are fermionic, so they must satisfy the anticommu- 
tation relations below [2T] . 

{Qa,Qi3} = {ga,g^} = o, (2.3) 

{Qa,Q,} = 2a^P^, (2.4) 

[Qa,Pf.] = [Qp,P,]=0, (2.5) 

[Q^,M^,] = i{aniQfi- (2.6) 

Here the SUSY generators are 2-component (Weyl) spinors, their indices and d,/? 
equal 1 or 2, and M^^ are the four-momentum and angular momentum operators, and 
cr^ and a^^ are defined in appendix |Al We have used the simplest case which is called 
= 1 SUSY, where A^ denotes the number of SUSY generators. 
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2.2.2 Superfields 

In order to describe the SUSY transformations, we introduce the superspace [20] which is 
different from the normal space by adding two new "Grassmann" coordinates 6a{a = 1, 2) 
and 90{P = 1,2). They are anticommuting 

{0a,efs} = {9^,9p} = {9^,9^}=0. (2.7) 

A "finite" SUSY transformation in the superspace {x^,9,9) can be defined as 

G{x,9,9) = exp[i{9Q + Q9-Xf,P^')]. (2.8) 

Here the indices have been dropped. The superfields $ should be functions of 9, 9 and 
such that they transform under the SUSY transformations as follows, 

G{y,ri,f])^{x,9,9) = ^{x + y - ir]a9 + i9af],9 + r],9 + f]) , (2.9) 

where rj and f] are again " Grassmann" variables. Considering infinitesimal SUSY transfor- 
mations on the superfield 

6s{v,m{x,G,G) = t{r]Q + Qn-x^P'')^{x,9,9) 

+ r/^ - 2 {w,0 - 9a,v) d'^] <!>{x, 9, 9) , (2. 10) 

one obtains the representation of the SUSY generators, 

= ^(J^-^/^^) ' (2-11) 



Qa = ^{4r + '(^^<J^9j . (2.12) 



The SUSY covariant derivatives anticommute with the SUSY transformations (12.10p . i.e. 

D^{6s<l>) = -SsiD^<^) . (2.13) 
This gives the expressions for the SUSY covariant derivatives. 

Da = + (2-14) 

Because the SM fermions are chiral, the first superfield we need is a chiral superfield, 
which can describe the left- and right-handed fermions as well as their superpartners. The 
left- and right-handed chiral superfields are defined via 

D^L = 0, (2.16) 
D^R = 0, (2.17) 
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respectively. The definition of tlie SUSY transformation, SUSY generators and SUSY 
covariant derivatives in the L(R) representation are similar to the definitions above [3]. 
The L(R) representation of the SUSY covariant derivatives are 



d „.„ _ ^ d 



Dl = — + 2ta'^9d^, Dl = ^, (2.18) 



Dr = -^-2iea''da, Dr = —. (2.19) 

do de ^ ' 

From fl2.16p we know that is independent of ^, therefore we can expand it as, 

$i(x,^) = <P{x) + V2d^ij^{x)+d^e^e^pF{x), (2.20) 

where Saf^ is the anti-symmetric tensor in two dimensions. From (12.71) we know that = 0, 
therefore the terms with three or more factors of 9 vanish. The field is a complex scalar 
field, and ip is a. left-handed Weyl spinor. They are called superpartners to each other. 
The field F is an auxiliary field, which we can get rid of by using the equations of motion. 
Under SUSY transformations the "component" fields (j),ip,F transform as 

6s{v,v)(j) = %/2r/7/>, (2.21) 
Ssiv,v)iJ = V2F + iV2a^f]d^(t), (2.22) 
5s{v.v)F = -tV2d^tPa^T] , (2.23) 

which shows that SUSY transforms fermions into bosons and the variation of the F field 
is a total derivative. 

The chiral superfield contains the spin-0 bosons and spin-| fermions. One also needs 
superfields which describe the spin-1 gauge bosons of the SM. Hence we introduce a vector 
superfield V{x, 9, 9) which satisfies 

V{x,9,9) = V^{x,9,9). (2.24) 
In component form it can be written as, 

V{x,9,9) = C{x) +t9x{x) -i9x{x) 

+ ^99[M{x) + iN{x)] - ^99[M{x) - iN{x)] 

- 9a^9v^{x) + i999[X{x) + ^^^d^x{x)] 

- im>^ + IcT^^d^xix)] + ^9999[D{x) + ^d^d^^Cix)] , (2.25) 

where C, M, N, D are real spin-0 fields, Xy ^ ^ire Weyl spinors and is a spin-1 gauge 
field. Under the non-Abelian supersymmetric gauge transformation the vector superfield 
transforms as 

e^gV _^ ^-^2gAt^2gV^^2gA^ ^2.26) 
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where A(x, 9, 9) is a chiral superfield, g is the gauge couphng, and V = V""T", A = A°T" 
where T° are the generators of the non-Abehan gauge group. For the Abehan case, the 
transformation simphfies to 

V — > V + i{A-A^). (2.27) 

One can choose the Wess-Zumino gauge [22]. where 

C{x) = x{x) = M{x) = N{x) = , (2.28) 

leaving the gauge field f^, its superpartner A (gaugino) and the auxiliary field D. In this 
gauge, the vector superfield is expressed as 

V{x, 9, 9) = -9a^9v^{x) + i999X{x) - i999X{x) + ^9999D{x) . (2.29) 

Similarly to the chiral superfield, one finds that the D component of the vector field trans- 
forms into a total derivative under SUSY transformations. This is very important for the 
construction of the SUSY Lagrangian. 



2.2.3 Construction of the SUSY Lagrangian 

In this section SUSY Lagrangian is constructed in the notation of superfields. The principle 
of construction is that the action should be invariant under SUSY transformations, i.e. 
6s J d^xC{x) = 0. It is satisfied when C transforms into a total derivative under SUSY 
transformations. From the discussions above we know that the highest components of 
chiral (F-term) and vector (D-term) superfields satisfy this demand. They can be obtained 
with integration over the Grassmann variables 9 and 9 via the definition. 



2.30) 



J d9a — , J 9ad9f3 — 5af3 ■ 

Therefore, the general form of the SUSY invariant Lagrangian is, 

C = I d^9CF + [ d^9d^9CD , (2.31 



where the Lagrangian densities Cp and Cjj are chiral and vector superfields, respectively. 

In order to obtain the explicit expressions for the SUSY invariant Lagrangian, we 
consider the product of left-chiral superfields. From the definition of the chiral superfield 
in fl2.16p . it is easy to see that the product of left-handed superfields is always a left-handed 
superfield. Therefore it can take the role of Cf- The superpotential iy(<l') (corresponding 
to Cp) is introduced as 

iy($i) = ^Ai<l>i + i^m,^$i$j + (2.32) 

i ij ijk 
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where $i are left-chiral superfields, and rriij and yijk are totally symmetric constants. 
Because of the renormalizability constraint the superpotential does not contain terms with 
four or more powers of superfields. The F-component of the superpotential can be written 
as 

i ij ' ijk 

This gives a fermion mass term and Yukawa interactions. 

Considering the product of the left-chiral superfield and its conjugate, one finds that 
is a vector superfield. Hence, it can take the role oi Lr, in (12.311) . The D-component 
of is obtained by 

y" dHdm^^^, = F,F* + d^<f),dy* + tiJ.a^d^iJ. , (2.34) 

which contains the kinetic terms for the scalar component (pi and the fermionic component 
ipi. Therefore the SUSY invariant Lagrangian can be written as 



d^dW{^i)+h.c. 



(2.35) 



In the discussions above, only the SUSY invariant Lagrangian is considered. But this 
is not sufficient. In order to describe the SUSY theory, we need a SUSY Lagrangian 
which is not only SUSY invariant, but also gauge invariant. Considering the non-Abelian 
supersymmetric gauge transformations, the vector superfield transforms as in fl2.26p . and 
the chiral superfield transforms as $ — 6"*^^"^$. The kinetic terms for the chiral superfield 
can be written as 

d^od'^e <i>+<i>, ^ / d'^ed^mte^^^^, 



g<t)*D(t), + F^ , (2.36) 

where the gauge-covariant derivative is defined by 

D, = d^ + igv;T'^ , (2.37) 

and are the generators of the gauge group. The first two terms of the second line in 
fl2.36p are the kinetic terms for scalars and fermions, and they also describe the interactions 
of matter fields (fermions and scalars) with the gauge fields. The couplings of fermions and 
scalars with gauginos can be found in the third term of the second line, while the third 
line contains the auxiliary fields D and F. 
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The kinetic term for the vector superfield is produced with tlie lielp of the field strength 
tensor, which is defined by 



W„ = ~DD[expi-2gV)D^expi2gV)] 



(2.38) 



where D and D are SUSY covariant derivatives and Wa is a left-chiral superfield because 
DoWa = 0. It is easy to show that the product WaW" is gauge invariant. Therefore its 
F-component can appear in the Lagrangian, 



(2.39) 



where 



F^u = dX - duv; + r'%y, (2.40) 

and /°*^ are the group structure constants. The kinetic terms for the gauge fields and 
gauginos fields, as well as the coupling of the gauginos to the gauge bosons are contained 
in ([239]). 

The Lagrangian invariant under SUSY and gauge transformations can be written as 



d'^9W{<!>i) + h.c. 



(2.41) 



The auxiliary fields can be integrated out with the help of their equations of motion. 



dC 

dC 



Da = -g^(t)*Ta(t)i. 



(2.42) 
(2.43) 



Substituting these expressions for the F- and D- fields into the SUSY Lagrangian (12.411) . 
the auxiliary fields in the Lagrangian disappear. The scalar potential can then be written 
as 



V 



'dW{(t)) 



+ 



9 



(2.44) 



where Vp and Vd are the F- and D-term of the scalar potential, respectively. 
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2.2.4 Spontaneous Breaking of SUSY 

The superpartners have the same masses as the corresponding SM particles, but no super- 
partners have been discovered yet. Hence SUSY must be broken. From the relation of the 
SUSY algebra (12. 4p one can derive 

H = P'^ = ^{QiQi + QiQi + Q2Q2 + Q2Q2) >0, (2.45) 

where H is the Hamiltonian. Suppose the vacuum state |0) is supersymmetric, which 
means 

Qa\0) = Qd|0) = 0, (2.46) 

then the vacuum expectation value is 

E^ac = {0\H\0) = 0. (2.47) 

This implies that SUSY is spontaneously broken ((5q|0) 7^ 0) if the vacuum expectation 
value is positive, i.e. E^ac > 0. It can be achieved if the scalar potential V does not vanish 
for any field configuration. Therefore we can break SUSY spontaneously with Da 7^ 
(D-term breaking) by the Fayet-Iliopoulos mechanism [23], or 7^ (F-term breaking) 
by the O'Raifeartaigh mechanism [21]. Unfortunately, none of these mechanisms works in 
the Minimal Supersymmetric Stand Model (MSSM). In the MSSM one introduces all the 
possible "soft-breaking" terms to the Lagrangian instead of assuming an explicit SUSY- 
breaking mechanism [25]. After these terms are introduced the quadratic divergences in 
6m'jj are canceled. The details of the soft SUSY-breaking will be explained in Section 

2.3 The MSSM 

2.3.1 The Field Content of the MSSM 

The minimal supersymmetric extension of the SM, which is called the Minimal Supersym- 
metric Standard Model (MSSM) [251 [27], keeps as few interactions and particles as possible. 
It means that the MSSM is also based on the gauge group SU{3)c x SU{2)l x U{1)y- It is 
not possible that one of the known SM particles is the superpartner of another one. Hence 
there must be a new superpartner for each known particle of the SM |27] . 

The SM chiral fermions and their superpartners are described by chiral superfields. The 
left-handed superfields can be arranged in S'[/(2)-doublets and the charged conjugates of 
the right-handed ones in S'?7(2)-singlets. For every generation we have five chiral super- 
fields: L, i?, Q, t/, and D. L and Q are the S'f/(2)-doublet chiral superfields which contain 
lepton(slepton) and quark (squark) doublets, respectively. E is the S'?7(2)-singlet chiral 
superfield containing lepton (slepton) singlets. There is no right-handed neutrino in the 
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SM, hence the corresponding superpartner does not exit. The S'L'^(2)-singlet chiral super- 
fields t) and D contain the up-quark (up-squark) and down-quark (down-squarks) singlets, 
respectively. In the SM gauge sector we introduce gluinos Aj(/c = 1, 2, ... 8) as the super- 
partners of the gluons G^, winos A"(a = 1,2,3) as the superpartners of the weak bosons 
ly^, and a bino as the superpartner of the hypercharge boson B^. They are described by 
the vector superfields V^, V"-, and v', respectively. 



Superfield 



Bosonic field 



Fermionic field 



SU,{3) SUl{2) Uy{1) 



L 

E 



sleptons 



leptons 



-1 
2 



Q 

u 

D 



Q 



squarks 
ul 



quarks 



di 



Q 



Ul 



u 



R 



di 



3* 
3* 



1/3 

-4/3 
2/3 



gluons 



gluinos A* 







ya 



weak bosons 



JL. 



wmos 



A" 







hypercharge boson 



bino 



A' 



Hi 

H2 



Hi 



H, 



Higgses 
Hi 



Hi 
Hi 
Hi 



Hi 



H, 



higgsinos 
Hi 



m 

Hi 
HI 



Table 2.1: Field Content of the MSSM 

The Higgs scalar boson has spin 0, hence we must describe it and its superpartner by 
a chiral superfield. But it turns out that one chiral superfield is not enough. One reason is 
related to the chiral (triangle) anomaly, which is proportional to Tr [F ^] , where Y denotes 
the weak hypercharge. In the SM, the chiral anomaly is canceled by the known quarks 
and leptons. However if one introduces only one chiral Higgs superfield, which contains 
chiral fermions with weak hypercharge y = lory = — 1, it leads to a nonvanishing 
contribution to the chiral anomaly, hence spoils gauge invariance. One therefore has to 
add a second Higgs doublet with opposite hypercharge. From Section 12.2.31 we know that 
the superpotential does not contain any products of left- and right-chiral superfields. Hence 
the Yukawa coupling term HQu^p^ cannot be replaced by HQu'j^ and the HQd'j^ and HLe'j^ 
terms cannot be replaced by HQd'j^ and HLe'j^, respectively. Therefore it is impossible 
to give masses to both up and down quarks if we only have one Higgs doublet. This is 
another reason why we need two Higgs doublets. Starting from the SM, we get the MSSM 
by introducing a superpartner to each known particles and adding another Higgs doublet 
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(with its superpartner) . The fields of the MSSM are summarized in Table [2?T] where the 
color and generation indices are suppressed. 



2.3.2 The MSSM Lagrangian 

As discussed in Sections 12.2.31 and I2.2.4[ the MSSM Lagrangian can be written as 

^ = ^susY + ^soft , (2.48) 

where Csusy is the SUSY-invariant term and Csoft is the soft SUSY-breaking term. All 
superfields of the MSSM were presented in Section 12.3. 11 The coupling constants of the 
gauge groups U{1), SU{2) and SU{?>) are denoted by g\ g and Qs, respectively. We define 
the field strengths of these groups as 

< = -\g'DDD^V', 

W» = -\DD[eM-'^9V)D^eM2gV)\, 

W^^ = -^DD[exp(-2(7,K)D,exp(2(7,K)] . (2.49) 

Here V = Yv', V = T^V, V, = T*^1/^ and Y, and are the generators of U{1), SU{2) 
and S't/(3), respectively. According to Section r2. 2. 31 the SUSY-invariant Lagrangian of the 
MSSM is written as 

Csusy = j dVO (le^''"+''^^ L + h9'v'+2gv^^ 

Q^g'V'+2gV+2gsVsQ _^ fj^g'V'+2gV+2gsVsfj _^ fj^g'V'+2gV+2gsVs jjj^ 
^^^g'V'+2gV+2gsVsjj^ _^ ^^^g'V'+2gV+2gsVs ^ 



^ ;W'W"' + -^W^W"- + -J-^WtWf + h.c. 



+ 



IQg'^ " 16^2 « lQg2 sa 

J d^Oe.j (^XuHiQ'U + \dH\QW + \^H{U E - ^H{Hi^ + h.c. 



,(2.50) 



where we have suppressed the color and generation indices. The last line in (12.501) is the 
F-component of the superpotential, where Sij is antisymmetric, i.e. ei2 = —S21 = 1 and fi 
is the mass parameter mixing the two Higgs fields. The Yukawa coupling constants A^, 
and Ae are 3x3 matrices in family space. 

The superpotential in (12.501) does not contain some terms which are gauge invariant 
and analytic in the chiral superfields. These terms are 

W = Eij [xUUE + X'UQW - ^'H\U^ + ybbt) , (2.51) 
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where A, A', A", and fi' are coupling constants. The first and second term in 12.511 violate 
the lepton (L) and baryon (B) numbers, respectively. B- and L-violating processes have 
never been seen experimentally, hence these terms are not included in the MSSM potential. 
In Section r2. 3. 3 1 we will analyze i?-parity which is related to B and L number conservation. 

In order to describe the MSSM completely, we need to specify the soft SUSY-breaking 
terms. Girardello and Grisaru [25] found out that the allowed terms are the scalar mass 
terms, gaugino mass terms and trilinear scalar interaction terms. The soft SUSY-breaking 
Lagrangian can be written as 

Csoft = -Ml(u{uL + d{dL)-M?^u^j,UR-Mj/jR 

- Ml (v{vL + e{eL) - MI/^cr 

- mfHlHi - mlHlH2 + {mleijH\Hi + h.c) 
+ ^ (Ml A' A' + MaA'^A'^ + MgA^A^ + h.c) 

- Sij (^XuAuHiQ% + XdAdH\Q^d\ + KA,H\U~e^^ + /i.e.) . (2.52) 

The first and second lines in fl2.52p consist of slepton and squark mass terms. The param- 
eters M?, M|^, Mj , M| and are 3x3 matrices in family space. In the third line 
we have the soft SUSY-breaking contributions to the Higgs potential with the Higgs mass 
parameters m\, m\ and m\. The gaugino mass terms are in the fourth line with the bino, 
wino, and gluino mass parameters Mi, M2, and M3. The last line in fl2.52p are trilinear 
scalar interaction terms with parameters A^, Ad, and Ag, which are 3x3 matrices in family 
space. Similarly to the superpotential in (12.501) . there are no terms which violate the B or 
L numbers in the soft-breaking Lagrangian. 

As discussed in Section 12.2.21 the superfields can be expanded in component form. The 
auxiliary fields (F- and D- fields) can be eliminated with the help of their equations of 
motion which are similar to (12.421) and (12.431) . Explicit expressions can e.g. be found in 
Ref.[2H]. 

For the quantization of the classical Lagrangian and for higher-order calculations we 
need to specify the gauge. A renormalizable 't Hooft gauge is chosen for this purpose. 
The principles can be found in Ref. [23 ED]- For the MSSM Lagrangian which contains 
the SUSY-invariant and soft SUSY-breaking terms, the gauge fixing Lagrangian can be 
written asl3T] 



- ^\d^W;: + iiwmwG^\' - {d^Glf (2.53) 

with arbitrary parameters ^a-, ^z, ^w, and ^g- Here A^, Z^, and W^^ are the physical 
electroweak gauge fields, are unphysical Higgs fields (Goldstone bosons). Since 
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Cfix involves the unphysical components of the gauge fields, one has to add the Fadeev- 
Popov ghost term Cgh for compensating these effects [30|. The complete Lagrangian of the 
MSSM is 

^ = ^SUSY + C,soft + C,fix + C,gh ■ (2.54) 

In the 't Hooft-Feynman gauge = = = = 1- The propagators for the 
unphysical Higgs fields then have poles at the masses of the corresponding physical particles 
and Z. 

2.3.3 i?-parity and Its Consequences 

In the SM the B and L numbers are conserved accidentally. In the MSSM, we can separate 
the most general gauge-invariant superpotential into two parts, one which conserves B 
and L numbers and the another one which violates them. Since B and L violation is 
not observed in nature, these terms must be suppressed or excluded. B and L number 
conservation together with spin conservation can be related to a new discrete symmetry, 
i?-parity^32j or equivalently matter-parity [33j. The -R-parity is defined as 

= (_l)3(B-L)+2.^ (2.55) 

where s is the spin of the particle. Thus, all SM particles have i?-parity Pr = 1, while all 
superpartners have i?-parity Pr = — 1. Usually i?-parity is assumed to be conserved in the 
MSSM, hence the superpotential W which violates i?-parity is forbidden. 

If -R-parity is assumed to be conserved, we have two important phenomenological con- 
sequences: 

• The lightest supersymmetric particle (LSP) must be stable. This means that super- 
symmetric particles other than the LSP must eventually decay into a state which 
contains an odd number of LSPs. 

• The supersymmetric particles can only be produced in even numbers in collider ex- 
periments. 

Moreover, a stable LSP has to be electrically neutral and interacts only weakly with ordi- 
nary matter, it can therefore make a good cold dark matter candidate [M] . 

2.4 The Physical Fields of the MSSM 

In this section the electroweak-symmetry breaking is discussed briefiy. The Higgs fields 
in the MSSM have eight real scalar degrees of freedom. When the electroweak symmetry 
SU (2) L X U{1)y is broken down to U {1)em, three of them turn into Goldstone bosons which 
are subsequently absorbed by and Z. The remaining five degrees of freedom form the 
physical Higgs bosons h^,H^,A^, and pS| [55] . After the electroweak symmetry is 
broken, fields with the same SU{3)c x U{1)em quantum numbers can mix. Therefore we 
have sfermions mixing and higgsinos and gauginos mixing into charginos and neutralinos. 
This will be studied in this section. 
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2.4.1 Higgs Bosons 

The scalar Higgs potential in the MSSM arises from the auxiliary F- and D-fields and the 
soft-SUSY breaking terms. It can be written as 

V = Vp + VD + Vsoft 

= {mj + \fi\^) hIh^ + {ml + HIH2 - mle.j {H{Hi + h.c) + 

\ + /) [hIh, - HlH,)\ ^^\HlH,\' . (2.56) 

We can take advantage of the arbitrary nature of the soft SUSY-breaking parameters 
ml, ml, and absorb into them, i.e. m^ + m\, ml + — >■ ml- Without losing 

generality, we can choose the vacuum expectation values of the Higgs fields as 

where Vi,V2 are non-negative. The electroweak symmetry SU{2)l x f/(l)y is broken down 
to U{1)em if these conditions are satisfied, 

ml + ml — 2m1 > , 

mjml - < . (2.58) 

The Higgs doublets Hi and H2 are decomposed in the following way: 

= I ^i + 75(</'? + ^X?) 



'1 



H2 = ( ^ 1 ^ • ON , (2.59) 

V ^2 + 7^(02 + %) / 

where 05 02 denote the neutral CP-even Higgs fields, Xi ^^id X2 denote the neutral 
CP-odd Higgs fields, and 0f and 0^ denote the charged Higgs fields. Inserting fl2.59p into 
the Higgs potential (12.561) . one obtains 



+ ^( X? xl)M,o [^h)+l{<^t <^t)M,. ( ) . (2.60) 

Here we only consider the linear and quadratic terms in the Higgs potential (12.561) . T^o 
and T^o denote the tadpoles of the Higgs fields 0? and 0°. They must vanish since the 
Higgs potential should have a minimum, 

-Amjvi + 4:mlv2 - {g^ + g'^) vi {vj - v^) 

2V2 

Amlvi - Amlv2 + (^^ + g'^) V2 (yl - v|) 
T,o = ^-^= = 0. (2.61) 
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In f l2.60p . M^o is the mass matrix of the neutral CP-even Higgs fields, 
M^o is the mass matrix of the neutral CP-odd Higgs fields, 

and M^± is the mass matrix of the charged Higgs fields, 

M ^ - ( mj + lig^vj + fvl) -ml-lg\^V2 \ . . 

"V -^l- ml + \ {g\l + ~g\l) ) ' ^^"^^^ 

with g"^ = g"^ + g''^, g^ = g"^ — g'^- These mass matrices can be diagonalized with the help 
of the unitary 2x2 matrices Ua, Up^, and Up^, respectively. 

These transformations result in the Higgs boson mass eigenstates: two neutral CP-even 
Higgs bosons and 

"° ^ -Uj^n, (2.66, 



2 / 

oTi^ rir^Uof^Tio K^or^Ti /CfO 



the neutral CP-odd Higgs boson and Goldstone boson G 



Upr. ( !h , (2.67) 



and the charged Higgs bosons 



G 



± 



'2 



In terms of the Higgs boson mass eigenstates, the scalar Higgs potential (12. 60p can be 
written as 

V = -n^h' - ThoH" + i ( ft" H" ) A/° (' ^° ) + (2.69) 
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Here Tyfi.Tno denote the tadpoles of the physical fields 



l'2 



From fl2.65bl) one gets, 



(cot /5 cos^ + tan /3 sin^ /?„ + sin 2/?„) , (2-71) 



where we have apphed fl2.61l) . Using the physical parameters e, Ow{,sw = sin 9w,cw = 
cos6'iy), tan/?, T^o, Tfjo, mz, m^o instead of the parameters vi,V2, g'^, g''^ , ml, ml, ml via 
f l2.7ip and the relations, 

tan/? = — , SH/ = l-m^/m|, e = g'cw = gsw, = hgl + gl){vl + vj) , (2.72) 

we obtain 

([/„M<^of/^)^^ = m|sin2(a + /5) + m^ocos2(a-/5)/cos2(/5-/5„) 

+ T^o cos(a - /5) sin^(a - /5„)/ cos^(/5 - 

2m^SH/Ciy 

e 1 

-T/jO- sin(a; — (cos(2q; — /5 — /3„) 



2m^svi/Ciy 2 

+3 cos(/3 - /5„)) / cos2(/5 - , (2.73a) 
([/aM(^of/J^)^2 = ""^1 sin(a + /?) cos(q; + /5) + m^o sin(Q; — /5) cos(a — (3)/ co^{(3 — (3n) 

+ 7^ sin(a - (3) sin^(a - /5„)/ cos^(/3 - /3„) 

2m^SH/Ciy 

-T/iO cos(a — (3) cos^(a — /?„)/ cos^(/3 — /3„) , (2.73b) 



2mzSwCw 

{U^M^oUl)^^ = mlcos\a + (3)+m\osm\a-f3)/cos\f3-f3n) + 

^ cos(a - Pn) (cos(2a -/?-/?„) 

-3cos(/?-/3„))/cos2(/?-/50 

- Tho sin(a - /3) cos^(a - /3„)/ cos^(/5 - /3„) , (2.73c) 

2mzSvi/Cvi/ 

(f/^„M^ot/^J = -mio tan(/3 - /3„) - ^^^^^^^ (Tj,o sin(a - /?„) + 

T/jO cos(a -/?„)) /cos(/5 - /3n) , (2.73d) 

f/^„M^of/^J^^ = m^otan2(/3 -/?„) + 

g 

{-Tho cos(a + /5 - 2/5„) 

2mzSvyCvy 

+Tho sin(a + /3 - 2/3„)) / cos2(/3 - /5„)) , (2.73e) 
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(Ui3,M^±Ul^ = m^o+m^, 

(Up^M^±Ul^^^ = - (m^o + m^) tan(/5 - 

g 

- 7^ {Tro sm{a - (5c) + T),o cos(q; - (5c)) / cos{(5 - (5c) , 

(UfiM^±Ul}j^^ = {m\o+ml,)iQ.Ti\(5-(5c) 

-Tho cos{a + (5- 2(5c)/ cos" {(5 - (5c) 



2mzSwCw 



+ Tho sm{a + P- 2(5c)/ cos^{(5 - (5c) . (2.73f) 

2mzSwCw 

The matrices UaM^oUl, Up^M^oUp^ and Up^M^±Up^ should be diagonal at tree level. This 
leads to the following conclusions, 

(5^^ (5^^ (5^ tan 2a = tan 2(5 ^^ ^ "^f , -- < a < J . (2.74) 

mj^ — m| 2 2 

Therefore, the tree-level masses of the physical Higgs bosons are 



^ {^\o + m| - ^^ {m\o + m|)2 - Am\om\ ^o^^ 2/3^ , (2.75a) 
^ (^m^o + m| + ^ (m^ + m|)2 - 4m^om| cos^ 2/3^ , (2.75b) 



m^± = m^o + . (2.75c) 
In the 't Hooft-Feynman gauge, the masses of the Goldstone bosons are 

777.(30 = 7772, = ■ (2.76) 
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2.4.2 W and Z Gauge Bosons 

In analogy to the SM, the physical gauge bosons PF^, and A^ are obtained from the 

electroweak interaction eigenstates PV^(a = 1,2,3) and via the following definition, 

A, = cwB^ + swWl, (2.77) 
and correspondingly their superpartners transform as 

= cw\' + sw>?- (2.78) 



The SU{2) X f/(l) covariant derivative is 



.Y 
2 

+ ^T+W^ + ^T-W- + 
^ [T^ - Qsl^] Z^ + leQA^ . (2.79) 

Cw 



D, = d. + tgT'^W^ + tg'-B, 



Here we introduced the new operators 

= T^±iT^, 
Q = + (2-80) 

where Q is the charge operator. Hence the kinetic terms for the Higgs fields are 

JO-maasMneUc = {D^^H,)^ {D,H,) + {D'^H^)^ {D.H^) . (2.81) 

After gauge symmetry breaking the Higgs fields acquire their vacuum expectation values 
as in fl2.57p . From fl2.8ip we can obtain the masses of the gauge bosons. 



9' 

= 2^ 



{vl + v^,) = l{9' + 9"){vl + vl), (2.82) 



and the masses of the gluons and photons are zero. This is consistent with the predictions 
of the SM. 
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2.4.3 Fermions 

Considering the Yukawa coupling terms in the superpotential in the SUSY Lagrangian 

jCvukawa = -Si, {XuHiQ'u'j, + Xd^Q^ d'j, + XeHlU 6%) + h.c. , (2.83) 

the fermion mass terms can be obtained when the Higgs fields get their vacuum expectation 
values, 

Cf-mass = -XuV2Ulu%- XdVidid^- XeVieLe% + h.c. . (2.84) 
Introducing four- component Dirac spinors, 

«=("0. e=('tM. (2.85) 



the fermion mass terms can be written as 

Cf-mass = -XuV2UU - XdVidd - XeViee . (2.86) 

Hence the fermion masses are 

rriu = XuV2 , rrid = XdVi , rrie = XeVi , (2.87) 

and the Yukawa coupling constants can be written as 

\ _ _ \ _ riidg X _ rrieg , , 

A^j — — — ^= , Ad — — — ^= , Ag — — — ^= . (^z.oo ) 

V2 v2sin/?mvi/ Vi \/2 cos (3ra\Y vi \/2 cos (3mw 

Here we have used tan l3 = ^ and = ^ ("^i + "^D- 

2.4.4 Sfermions 

The mass terms of the sfermions arise from the soft SUSY-breaking Lagrangian, the aux- 
iliary F- and D-fields. When the electroweak symmetry SU{2)l x f/(l)y is broken down 
to U{1)em, the Higgs fields get their vacuum expectation values, the sfermion mass terms 
can be written as 

^ioft = -Ml[u{uL + d{dL)-Mi^u^j,UR-Ml/jR 

- Ml (z>iz>L + ^gi) - Mle^CR 

- (^XuAuV2UlUr + XdAdVididR + Ae^ef le^e^ + h.c.^ , (2.89) 



^Aux-F = -Xlv\ (e[eL + - A>i (d[dL + d}j^dR^ 

- Xlvl {u\ul + U^rUr^ 

+ {XeV2^ie^ReL + XdV2fidRdL + XuVifiu'j^UL + h.c.^ , (2.90) 



24:CHAPTER 2. THE MINIMAL SUPERSYMMETRIC STANDARD MODEL (MSSM) 



'-Aux-D 



vl-vl 



-m^ cos 



25 E 



~fl{Tl-Qf,s'w)~fL+~fiQfnSyR 



(2.9i; 



Here we have applied ^ ^^^a 



m| COS 2/3. T^, F and Q are the weak isospin, weak 



hypercharge and electric charge of the fermions, respectively. They satisfy the relation 



fR 



Q = T^ + ^. For the right-handed fermions = 0, so Qf^ — ^ . 

Adding together the Lagrangian above, the sfermion mass terms of the Lagrangian can 
be written as 



f-mass - -~^{M,UL - {fl H) ( £ 



(2.92) 



After replacing the Yukawa coupling constants Ae, A^, and A^ with the expressions in (12.881) . 
we get the sfermion mass matrix M. p 



'm) + M?^ + m| cos 2/3(T| - QfS 
mf{Af — /ik) 



w) 



mf{Af-fiK) 
mj + M?^ + m| cos 2(3Q fS 



2 / 5 

W 



(2.93) 



where = M| , = M? , and the parameter k is defined as k = cot /3 for up- type 
squarks and k = tan (3 for down-type squarks and sleptons. The sneutrino mass matrix 
Mj> for a given flavor is 1-dimensional, and only has the left-handed entry of (12.931) . 



Ml + -m| cos 2(3 . 



(2.94) 



The sfermion mass matrix (12. 931) can be diagonalized by a unitary 2x2 matrix U ^ 

^2 



UjMjU} = Mf 



ji 

771% 
72 



(2.95) 



with the mass eigenvalues m% , m^- , 

h h 



m 



fl.2 



^ {Ml + M?J + mJ + ir|m| cos 2(3 

1 
2 



T-\/ [Ml - Ml + ml cos 2P{Tf - 2Q;.^)] + Am}{Af - f^n)^. 



(2.96) 
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The sfermion mass eigenstates are given by 



/l \ TT I f' 



f2 J '^ni 



where the matrix U? is determined by fl2.95p . 



IR 



Ufi f , (2.97) 



2.4.5 Charginos and Neutralinos 

The mass terms of the higgsinos and gauginos arise from the SUSY Lagrangian 

- V2 [^#1 (^gX + ^^'A'^ + i&2 [g\ + \g'\'^ + /^-c] , (2.98) 
the superpotential 

Si, (^fiHiHi + h.c^ , (2.99) 

and the soft-breaking terms 

^ (MiA'A' + MaA^A'^ + h.c.) . (2.100) 

The charginos are a mixture of the charged higgsinos and gauginos. We introduce 2x1 
matrices 

for the charged higgsinos and gauginos, where X± is defined via fl2.78p . Apply (12.1011) to 
the mass terms of the higgsinos and gauginos Lagrangian (12.981 12.99[ 12.1001) . one obtains 
the mass term of the charginos after the Higgs fields acquire their vacuum expectation 
values, 

C^^-mass = -[^Ij''^X^'^+^''^X^^''], (2.102) 

with the mass matrix 

^ ^ f M, V2m^ sin/? \ _ 3 

\ V 2 mw cos p II J 

The mass matrix (I2.103P can be diagonalized by two unitary 2x2 matrices U and V, 

UXV^ = M^,= i j, (2.104) 
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with the eigenvalues m 



1 1 

171-+, 171- + 



2\2 



+ 4m^ cos^ 2P + 4:ml^{Ml + /i^ + 2 /i M2 sin 2/3) 



Hence the chargino mass eigenstates can be written as 



(2.105) 



(2.106) 



with the unitary matrices U and V which are determined by fl2.104p . The four- component 
chargino spinors are defined as 



'X.i 
'X.i 



1, 2. 



(2.107) 



Similarly to the chargino case, the neutralinos arise from mixing of the neutral higgsinos 
and gauginos. In the following we introduce a 1 x 4 matrix 



( -tXA , -tXz , Hi , HI ) 



(2.108) 



to the Lagrangian of the higgsinos and gauginos (12.981 I2.99[ I2.100p . where Xa,Xz are 
defined in (12.781) . After the electroweak symmetry is broken, one obtains the mass term of 
the neutralino Lagrangian 



•^x"— mass 



^qT, , —0-1 



(2.109) 



with the mass matrix 



Y 



Ml —mz sw cos (3 mz sw sin P \ 

M2 mz cw cos 13 —mz cw sin [3 

—mz Sw cos [3 mz cw cos /5 —fi 

\ Sw sin P —mz cw sin /3 — /i J 



(2.110) 



In order to diagonalize the mass matrix and get the neutralino mass eigenstates, the fol- 
lowing transformations must be performed: 



/ m^o 



N*YN^ 



m^o 

X2 










V 



m.r.0 
X3 



\ 







(2.111) 



(2.112) 
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where is a unitary 4x4 matrix and x° are the neutrahno mass eigenstates. One of the 
eigenvalues in fl2.11ip might be negative if the matrix N is reaL Therefore, this matrix 
should be complex for positive neutralino masses even if all the elements in Y are real. 
The neutralinos are Majorana fermions, their Majorana spinors are defined by 

X° = (l^), ^ = 1,...,4. (2.113) 

We have four neutralinos Xi, X25 Xs? xl- They are labeled in ascending order, m^o < m^o < 
m^o < m^o. The lightest neutralino Xi is the only MSSM particle which can make a good 
cold dark matter candidate. 



2.4.6 Gluinos 

The gluinos are the superpartners to the gluons G^. Their mass terms arise from the 
soft-breaking terms of the Lagrangian, 

£ = i + /i.e.) . (2.114) 

Since the gluinos are Majorana fermions, their Majorana spinors are defined by 

~9s = (I) (2-115) 

with the Majorano mass M3. 



2SCHAPTER 2. THE MINIMAL SUPERSYMMETRIC STANDARD MODEL (MSSM) 



Chapter 3 

Renormalization of the MSSM 



In order to calculate higher-order corrections, one must renormalize the parameters and 
fields of the MSSM. Several approaches for the renormalization of the MSSM have been 
developed [HI [ISl EZl ESI ESI SO] • Here we follow the strategy of Refs. [HI [13] , i.e. we employ 
on-shell renormalization. We assume here that all relevant parameters are real quantities. 
This amounts to the assumption that the soft supersymmetry breaking terms conserve 
CP. In this chapter the basic ideas of regularization and renormalization are presented 
briefly. The DR and on-shell renormalization scheme are introduced and explicit on-shell 
renormalization conditions for the different flelds are formulated. The renormalization of 
the SM sector, the chargino and neutralino sector, the sfermion sector, and the Higgs sector 
are discussed in detail. 



3.1 Principles of Regularization and Renormalization 

At tree level, the parameters of the Lagrangian are directly related to the physical quan- 
tities. In higher-order perturbation theory, these direct relations are destroyed. Further- 
more, the divergent loop integrals make the calculation ambiguous. The theory has to 
be regularized for mathematical consistence. One usually employs dimensional regular- 
ization [?T] for higher-order calculations of the SM, where the integrals are calculated in 
D dimensions. Dimensional regularization preserves Lorentz and gauge invariance of the 
theory, but breaks chiral symmetry and SUSY. In supersymmetric theories, if we treat 
the vector fields in D dimensions, this will lead to a mismatch between the fermionic and 
bosonic degrees of freedom, hence SUSY is broken. In order to avoid the disadvantage of 
dimensional regularization, dimensional reduction was developed in Ref.[12], where only 
the momenta are calculated in D dimensions, while the fields and the Dirac algebra are 
kept 4-dimensional. It is commonly used in supersymmetric theories. 

After regularization the parameters in the original Lagrangian (the so-called bare 
parameters), which are physically meaningless, are different from the corresponding physi- 
cal quantities by UV-divergent contributions. These divergences cancel in relations between 
physical quantities. Therefore, one may replace the bare parameters by the renormalized 
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ones and the renormalization constants (counterterms) in the Lagrangian, 

9o — ^g + Sg, (3.1) 

where the renormahzed parameters g are UV finite and measurable, while the counterterms 
Sg are UV divergent and absorb the divergent parts of the loop integrals. They are fixed 
by renormalization conditions, which determine the relation between the renormahzed 
parameters and the physical quantities. One must choose a set of independent parameters 
in order to make predictions from the theory. 

Parameter renormalization is sufficient to get finite S-matrix elements when vave func- 
tion renormalization for external on-shell particles is included, but the off-shell Green func- 
tions are not finite. In order to obtain finite propagators and vertices, we must renormalize 
the field by the multiplicative renormalization 

V'o — > -y/SZ^ip i/j + ISZ^tjj . (3.2) 

After the parameters and fields are renormahzed, the Lagrangian can be written as 

^o(^o,^o) = C{g,ij) + SC{g,5g,ip,5'i/j), (3.3) 

where the renormahzed Lagrangian C has the same form as the bare Lagrangian Cq but 
depends on the renormahzed parameters and fields and SC contains the counterterms. 



3.2 Renormalization Scheme 

The renormalization constants can be chosen arbitrarily. Their divergent parts are deter- 
mined by the structure of the relevant loop integrals, while their finite parts depend on the 
choice of the renormalization conditions. Here we focus on two different renormalization 
schemes: 

• On-shell renormalization scheme 

The counterterms are fixed in a way such that the finite renormahzed parameters are 
equal to some physical quantities. 

• DR renormalization scheme 

The counterterms are defined such that they only contain the UV-divergent parts of 
the bare parameters or fields. The UV-divergent parts are proportional to A, which 
can be written as 

A = ^A__^g + iog47r, (3.4) 
where 7e is Euler's constant. 
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3.3 On-shell Renormalization Scheme 

The basic idea of the on-shell renormalization scheme is: 

• The counterterms for parameters of the physical particles are chosen in a way such 
that the renormalized masses are equal to the physical masses, which are the real 
parts of the poles of the corresponding renormalized propagators. 

• The diagonal entries of the field renormalization matrix are fixed by the requirement 
that the real parts of the renormalized propagators have unity residues. 

• The renormalized IPI two-point function (the inverse of the renormalized propagator) 
is diagonal for on-shell external particles. This determines the non-diagonal entries 
of the field renormalization matrix. 

In the on-shell renormalization scheme all renormalization conditions are formulated for 
on-shell external particles. Now we come to the on-shell renormalization conditions for 
different types of fields. 

3.3.1 On-shell Renormalization Conditions for Scalars 

Consider the scalar one-particle irreducible(lPI) diagram iTj{p^). Via the Dyson summa- 
tion the full propagator can be written as 



iA(p^) ^ — + — O" + — Q-Q" + 

pz — 777,^ — — — — 



p2 _ _|_ Xl(p2) 



(3.5) 



If there are n mass eigenstates 0j(i = 1, • • - n), the renormalized IPI two-point function 
can be written as 

^ij{p) = i{p'^ -mf)5ij + i%j{p^), (3.6) 

Eij(p^) is the renormalized IPI self-energy. 

The on-shell renormalization conditions require that the poles in the renormalized prop- 
agators occur at p^ = m| and the renormalized IPI two-point function is diagonal, which 
are equivalent to 

Ref,,(p)|^,^^, =0 ^ Re%^{m]) = Q. (3.7) 
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Re takes the real parts of the loop integrals in the self-energies. It can be replaced by 
Re if all the relevant couplings are real. Furthermore the real parts of the renormalized 
propagators have unity residues, 

lim ^ ^ Retiiip) = i =^ ReS^(m2)=0, (3.8) 
where E^(mf) = ^t^ip'] 

3.3.2 On-shell Renormalization Conditions for Fermions 

Similarly to the scalar case, the renormalized IPI two-point functions for the fermions can 
be written as 

ff^.(p) = z5,,(^-m,)+zS,,(p2). (3.9) 
Sjj(p^) is the renormalized IPI self-energy for the fermions, which can be decomposed via 
t,{p') = ^..^Sj.(p2)+^.;^Sg(p2)+.;^Ef (p2)+^«Sf (p2). (3.10) 

where ul,r = (1 =F 75)/2. 

The on-shell renormalization conditions for fermions are 

ReWp)fi,{p)L, = 0, /x.(p')Reff,(p')l„,2_. = , (3.11) 



lim Ret{-{p)fii{p) = ifii{p), lim fii{p')Retl{p') = ^f^iip')^ (3-12) 

here fi{p),fi{p) are Dirac spinors of the external fermion fields. The diagonal equations of 
(13. lip ensure that the renormalized fermion masses are the poles of the corresponding prop- 
agators, while its non-diagonal equations make the renormalized IPI two-point function 
diagonal for on-shell external particles. The conditions that the renormalized propagators 
have unity residues are satisfied by (I3.12p . These on-shell renormalization conditions are 
translated into the relations between the renormalized self-energies for the fermions, 

mjReEj;.(m2) +ReSf/(m|) = 0, (3.13a) 
mjReSg(m2) +ReSg^(m2) = 0, (3.13b) 

ReE,^(m?) + ( Ret^^' (mj) + ReSf (m? 



+m, (ReSf (m2) + ReS,f (m^)) = 0, (3.13c) 
R~eE,^(m^) + (ReE^^VD + ReS.f (m^)) 

+mi(Re±fi^'{m^i) + Ref:fi^'{m^,)) = 0. (3.13d) 
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3.3.3 On-shell Renormalization Conditions for Gauge Bosons 

Using the 't Hooft-Feynman gauge, the renormahzed IPI two-point function for the gauge 
fields can be written as 

r%ip) = -^9,Ap'-^^)^^^+^[9,.'^)mp')-^^-f'mp')■ (3-14) 

Here E*^ (p^) is the renormahzed self-energy for the gauge fields, and the indices T and L 
denote the transversal and longitudinal parts, respectively. The on-shell renormalization 
conditions require that the renormahzed masses are equal to the physical masses and the 
renormahzed IPI two-point function is diagonal when the external particles are on their 
mass shell, 

Ret%{p)e'<{p) |^,_, = , Ret%{p)e'^{p) |^,_, = , (3.15) 
and the renormahzed propagators have unity residues, 

lim^ -^^Ref;V(p)er(p) = -^e,Jp) , (3.16) 

where are the polarization vectors and satisfy p^e^{p) = when = m?. Applying 
f l3.14p to the on-shell conditions 03.15^ 13.161) . one obtains the on-shell conditions for the 
renormahzed self-energies, 

ReE^ (m2) = , ReE^'(m^) = . (3.17) 

Note that the longitudinal part of the gauge boson self-energies are dropped since they are 
always finite. 

3.4 Renormalization of the SM-like Sector 

In the SM the input parameters are chosen to be the electric charge e, the fermion masses 
ruf, and the masses of the W and Z gauge bosons. The on-shell renormalization of the 
Standard Model has been performed in [TT], H3]. Here we follow the conventions of Ref. [11] . 

3.4.1 Fermion Sector Renormalization 

The bilinear part of the Lagrangian for SM fermions is 



(3.18) 
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where / = iy,e,u,d denotes the four- component Dirac spinors and i = 1,2,3 denotes the 
generation index. The fermion masses and fields are renormahzed via 

rrif + 6mf , (3.19) 



ruf 



6ii + 



f,R 



3 1 



(3.20) 



where 5m j is the counterterm for the fermion masses mj, and 5Z(-^ and 5Z(-^ are the field 
renormalization constants for the left- and right-handed fermion fields, respectively. Ap- 
plying the transformations (13.191 [3l20l) to (13.181) . one obtains the counterterm Lagrangian, 

5ZIi^ul + SZffuR 



5C 



/, - [mj^ (dZl:"^ + + 6m f] /,/, . (3.21) 



In general the renormahzed self-energy is equal to the unrenormalized self-energy Ejj(p^) 
plus the corresponding counterterms, which are the derivatives of the counterterm La- 
grangian 5C with respect to the fields fi and /j. 



Therefore the renormahzed self-energies for the SM fermions can be written as 



(3.22) 

(3.23a) 
(3.23b) 



Mfij^) + ruf^ {5Z(f + 5Z(f) + Smj^ . (3.23c) 



Assuming the CKM matrix as an identity matrix Vij = Sij, all the field renormalization 
constants and the self-energies are diagonal, 

, for i ^ j . (3.24) 

Applying the renormahzed fermion self-energies (13.231) to the on-shell renormalization con- 
ditions (13.131) . one can fix the renormalization constants 



szl:' 



^ (ReEr(mJJ +ReS^^(m|)) 



2 Y It y-'h 



m 



(m^J-m|(ReS^'^'(mJJ+ReS 



f,R' I 2 ^ 



2mj-ReS^''^ (m 



2 1 



5z!f = -ReS^'^(mJJ - m]^ (ReSf (mjj + ReS^'«'(mJj) 



2m f- ReS 



(3.25a) 

(3.25b) 
(3.25c) 
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3.4.2 Gauge Sector Renormalization 

The bilinear part of the Lagrangian describing the gauge fields is 

'(^ = W~{p'-ml,)W^>' + Z^{p'-ml)Z^^ + A,p^A>^. (3.26) 

The gauge sector is renormalized via the transformations 

ml^r — > + (5m^ , (3.27a) 

m| — > m| + (5m| , (3.27b) 

— > + ^6ZwW^ (3.27c) 



Z \ / 1 + l5Zzz \5ZzA \ ( Z 



A J V '2^Zaz 1 + kSZAA J V ^ 



(3.27d) 



After the renormalization one obtains the Lagrangian which gives us counterterms 

6C = W-[6Zw{p'-ml,)-6mlr]W+^ 

+ [6Zzzip' - ml) - 5ml] Z^ + A^A^SZaa 

+ ^A, [6Zza{p' - ml) + 6Zazp'] . (3.28) 
The renormalized self-energies for gauge fields can be written as follows, 

t^ip') = E^{p') + 6Zw{p'-mly)-6mlr, (3.29a) 

E|^(p2) = S|^(/)+5Zzz(p'-m|)-5m|, (3.29b) 

t^^ip') = T.^\p^)+5Zaav\ (3.29c) 

t^'{p^) = Y.f{p') + ^ {5Zza{p' - ml) + 6Zazp') ■ (3.29d) 

Using the on-shell renormalization conditions (13.171) . we can determine the renormalization 
constants for the gauge sector of the SM, 

Ret^{mly) = ^ = ReJ:^{ml,) , (3.30a) 

ReS|^(m|) = ^ 5m| = ReE|^(m|) , (3.30b) 

Re±f{ml) = =^ SZaz = -'^^^^^^^ , (3.30c) 

mz 

ReE^^(0) = =^ 6ZzA = ^^^^f , (3.30d) 

ReE5^'(m^) = ^ 6Zw = -ReE^'(m^) , (3.30e) 

ReS|^'(m|) = =^ 5^^^ = -ReE|^'(m|) , (3.30f) 

ReS^^'(O) = ^ 6Zaa = -ReE^^'(O) . (3.30g) 
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The on-shell definition of the weak mixing angle Ow{sw = sm9w, cw = cos 9w) is [33] 



,2 _ 1 ^VK 

''w ~ ^ 2 • 
mi 



(3.3i; 



"Z 

Hence its counterterm is directly related to the counterterms of the gauge boson masses, 

6sw f 



sw 2 V ™i 

5cw _ If 
Cw 2 V m| 



(3.32) 



3.4.3 Electric Charge Renormalization 

The three-point function 7// vertex at one-loop level can be depicted as 



/, P' 




/, p 



-ieQf-f^ 



+ 



+ 




+ 




(3.33) 



The on-shell renormalization condition for the electric charge requires that all corrections 
to the 7// vertex should vanish for the on-shell external particles in the Thomson limit 
{p = p'), 



-ieQfu{p)-f^u{p) 



u{p)Tl^^{p,p)u{p) 

^ u{p) {PC<Jf{p.p) + 5Alff) u{p) = . 

The electric charge is renormalized via 

e — ^ (l + 5Ze)e. 



(3.34) 



(3.35) 



Together with the fermion field fl3.20p and photon field transformation (13.27dp one obtains 
the counterterm for the 7// vertex, 



+ 7t^iVf - af-f5)-6ZzA: 



(3.36) 
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where Vf and a/ are the vector and axial vector couphng of the Z boson to the fermion /. 
Inserting the Ward-identity 



d 



u{p)Alff{p,p)u{p) = -Qfu{p)^—J:ff{p)^u{p) -2afu{p)-f^ujLu{p) ,(3.37) 



rriz 



where 



u 



(p) \l-j:jf(p)] u{p) = -u{p)j, {SZ^'^^UL + SZf'^cun) u{p) , (3.38) 



dp' 



into the charge renormahzation condition (13.341) . we can fix the charge renormahzation 
constant dZ^.. Its explicit expression is 

5Ze = -\bZAA-\—^'Z'ZA, (3.39) 

where we have used the relation Vf — af = —^Qf- 

The fermion-loop contributions to the photon field renormahzation constant SZaa in 
(I3.39P give rise to large logarithm Inmj (/ denotes the light fermion). The fine structure 

2 

constant a = j;^ and the Fermi constant have the relation as 

n 1 + Ar ^ ^ 

where Ar summarizes all the radiative corrections to the muon decay [Tl| 1^ HHl HT] . One 
therefore can parameterize the Born matrix element by G^, i.e. 

2m^s^ I ^) ^ • (3.41) 

Combine the charge renormahzation constant (I3.39P and the Ar contributions in (13.411) . 
one obtains 

6Ze = 6Ze Ar 

2 

lc^^/m| Sml^^ l S^(Q)-^m^ 1 Sj^^(O) 1 

2s^ mi. 2 swCw mi. 2 

where Arbox denotes the box corrections to the muon decay. Thus the large logarithm 
Inrrif disappears in (I3.42p . 
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3.5 Renormalization of the Chargino and Neutralino 
Sector 

Renormalization of the Chargino Sector 

The kinetic and mass terms of the chargino Lagrangian in terms of the four-component 
chargino spinors xt "^^^ t)e written as 

^ = ^[#.. - (M^^h - iM^+)Jj] xt , (3.43) 

where MP+ = UXV~^ is the diagonahzed mass matrix of the charginos. In order to renor- 
mahze the chargino sector, we introduce the counterterm for the chargino mass matrix 

X — > X + 6X, (3.44) 



6X 



6M2 V2 6 {mw sin P) 

\/2S(m\Y cos/?) Sjj 



(3.45) 



in which the counterterms 6M2 and 6fi are determined in the chargino sector renormal- 
ization. The counterterm for the W boson mass has been determined in fl3.30al) and the 
renormahzation of tan /? will be discussed in Section 13. 7[ 

The chargino fields are renormalized via the transformations 



(3.46) 



where the field renormalization constants 6Z^, 6Z^ are general 2 x 2-matrices. Applying 
the transformations (I3.46P and (I3.44p to the Lagrangian (13.431) . one gets the counterterm 
Lagrangian 



-X- 



2 



(3.47) 



[{U5XV^y + t^M^, ' 5Z^* + Uz^'M^,\u:n 



According to (13.101) and (13.221) . the renormalized self-energies Sjj(p) for the charginos can 



3.5. RENORMALIZATION OF THE CHARGINO AND NEUTRALINO SECTOR 39 



be written as 



J:Up') + 1{6Z^* + 6Z^^ 



U6XV^ + -6Z^'^MP+ + -MP+6Z^ 
2 ^ 2 ^ 



(3.48) 



The counterterms 5M2 and 5^ are determined by renormalizing the two charginos via 
the on-shell renormahzation scheme. According to the on-shell renormahzation conditions 
for the fermions in Section I3.3.2[ we obtain the on-shell renormahzation conditions for the 
charginos, 



m- 



,Ret^Am^.+ )+Ret 



Ai Ai \ Aj Aj 

+2m^+i?eEf' (m^+) 



0. 



(3.49a) 
(3.49b) 

(3.49c) 

(3.49d) 



f l3.49cp and fl3.49dp . which make the renormalized chargino propagators have the residues 
1, fix the diagonal entries of the chargino field renormahzation matrices, 



5Zi 



SZH 



2m-+ReEf/'\m^.+ 



(i?eS^'(m^0 + i?eS,f(m2 

Ai \ Ai A^ 

I?) 

2 



2m,-,+ i?eEf^'lm 



m 



xi 



(3.50) 



(3.51) 



The diagonal equations of fl3.49ap and fl3.49bp . which ensure that the renormalized 
chargino masses are the poles of the corresponding propagators, determine the counterterms 
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5M2 and 

6M2 



S/i — 



- Ui2Vi2(m~+ [i?eE^2(^+) + Re^g{ml+)] + 2ReJ:li{ml+)^ 
+ 2 (C/12C/21 - C/iiC/22)T42l^22 5(V2mw sin f3) 

+ 2 t/i2t^22(^i2V^2i - VL1V22) diV2mw cos/3)] /A , (3.52) 
UnVn (m-+ [i^eE^^^l^j) + i?eE|(m2 + )] + 2i?eEf2^(m2 + )) 

- U21V21 (m-+ [i?eEf'i(m^+) + i?eEfi(m^+)] + 2i?eEff'(m2 + )) 
+ 2 f/nf/21 (V^i2V^2i - V^ii^22) ^(y2m^y sin/5) 

+ 2 (?7i2^72i - ^/ii^/22) ^11^21 S{V2mw cos /3)1 /A , (3.53) 



with A = 2([/ii[/22ViiV22 — ^12^21^12^1)- In contrast their non-diagonal equations, which 
make the renormahzed chargino IPI two-point functions diagonal for on-shell external 
particles, determine the non-diagonal entries of the chargino field renormalization matrices. 



szh = 



m. 



mU ReEUmU) + m~+ m-+ ReY.fAm'+) + m-+ i?eEg^(mL) 



+ ReE?^{mi+) - {USXV ' )ij - {U5XV 



m. 



mi: 
xl 



ml+ ReT.fAml+) + m-+ m~+ ReT.^^{m\+) + m-+ ReT^i^iml 



{U5XV^)ij - {USXV^)j^^ . (3.54) 



3.5. RENORMALIZATION OF THE CHARGINO AND NEUTRALINO SECTOR 41 



Renormalization of the Neutralino Sector 

In terms of the neutralino Majorana spinors the Lagrangian which describes the kinematic 
and mass terms of the neutrahnos can be written as 



1 ^ 



(3.55) 



In analogy to the chargino case, we introduce the counterterm for the neutralino mass 
matrix Y and the field renormalization constants for the neutralino fields by the transfor- 
mations 



Y 



Y + SY. 



^LXj 



^^RX°i' 



(3.56) 



(3.57) 



where the field renormalization constant 6Z^ is a general complex 4 x 4-matrix. One does 
not need to renormalize the left and right components of neutralinos independently due to 
the definition of their Majorana spinors in fl2.113p . The elements of the matrix 6Y are the 
counterterms for the parameters in the mass matrix (12.1 lUp . 



6Y 



( SMi — S {mz sw cos f3) 5 {niz sw sin (3) \ 

6M2 6 {mz cw cos /?) —6 {mz cw sin P) 

—6 {mz Sw cos P) 5 {mz cw cos P) —6fi 

\ 5 {mzsw sin P) — 5 {mzcw sin P) —djji / 



.(3.58) 



Applying the transformations (I3.56[ 13.570 to the Lagrangian (13. 550 . one arrives at the 
counterterm Lagrangian 



CT 



2 



, 6Z^^ MP + MP6Z^' 



( , ^ MP6Z^* + 6Z^^MP 
+ NSY^N^ + ^ 



(3.59) 
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here is the diagonahzed neutrahno mass matrix which has been defined in (12.1111) . 
Similarly to the chargino case, the renormalized neutalino self-energies can be written as 



X 7O I j\/fD _|_ j\/fD r 7O 

tfHp^) = EfHp') - I N*6YN^ + ^» ^ ^° I . (3.60) 



Obviously they obey the relations 



^Up') = ^Up') ' (/) = Ef (p^) , (p^) = Sj-(/) , (/) = (p^)t . 

(3.61) 

Only the counterterm for the parameter Mi in the chargino/neutralino sector is not 
determined so far. We can fix it by renormalizing one of the four neutralinos via the 
on-shell renormalization scheme. Conventionally, the lightest neutralino Xi is chosen for 
this task. From the previous discussions, the on-shell renormalization conditions for the 
neutralino sector can be expressed as follows. 



m^oi?eSg(m|o) + i?eSf/(m|o) = 0, (3.62a) 

m^oi?eSg(m|o) + i?eSf/(m|o) = 0, (3.62b) 
for {t ^ j) y{t = j = l) 
Ret^i{mlo) + 2mloRet^\mlo)+2m^oRetf/''{mlo) = 0. (3.62c) 

The diagonal equations (i = j = 1) of (I3.62ap and fl3.62bl) . which ensure that the renor- 
malized lightest neutralino mass is the pole of the corresponding propagator, determine 
the counterterm 6 Mi. In contrast their non-diagonal equations, which fix the non-diagonal 
entries of the neutralino field renormalization matrix, make the renormalized neutralino 
IPI two-point function diagonal when the external particles are on their mass shell. The 
diagonal entries of the neutralino field renormalization matrix are determined by fl3.62cl) . 
which make the renormalized neutralino propagators have the residues 1. 

Inserting the renormalized neutralino self-energies fl3.60p into the on-shell renormaliza- 
tion conditions, one obtains the expressions for the renormalization constants 6M1, SZ^^ 
and 6Z^, which are as follows. 



SMi 



i^ii 



2Nii [A''i3 5{'mz sin 9w cos (3) — Nn 5{mz sin 9w sin /?)] — 

2Ni2 [Nrs 5{niz cos Ow cos (3) - 5(m^ cos dw sin (3)] - N^^ 6M2 

+2Ni3Ni^ 6jj + m^oi?eSfi(m|o) + i?eSff (m|o)j , (3.63a) 
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= -i?eS,^(m|o)-2m^o[m^oi?eS,^'(m|,) + i?eSf'(m|o^^ (3.63b) 
2 ~ 

Xu -^u 

+ m^oi?eSg^(m|o) - m^oXN5Y N'^)ij - m^o^{N*5YN^)ij\ . (3.63c) 



3.6 Renormalization of the Sfermion Sector 



The kinetic and mass terms of the sfermion Lagrangian are 

^ = (/! fl){v'-Mf)(^h'^+ul{p'-M,^)f>,, (3.64) 

where Mf = UfMfU\ is the diagonahzed sfermion mass matrix. At one-loop level, the 
counterterms for the sfermion mass matrices Mj and Mj^^ are introduced via 

Mj Mj + 5Mj , ^ Mi>, + 5Mp, . (3.65) 

The elements of the matrix 5Mj and ^M^jj are the counterterms for the parameters in the 
mass matrices (12. 93|) . (12.941) . respectively, 



6Mf 



6M, 



5 Mi 

ij + 5C 



mf6Aj 



/l2 



6M 



+ ^ ( -'m-'z cos 2/3 



(3.66a) 
(3.66b) 



where 



5Cj_^^ = 5{m} + mlcos2P{Tf-Qfsl.)) , 

SCj^^ = 5(mJ + m|cos2/3Q/s^) . (3.67) 

In order to get finite Green functions, the field renormalization constants are introduced 
via the transformations 

''l + l-SZf] ( V (3.68a) 



1 

2 



^ [1 + Uz,^ ) i)i . (3.68b) 
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The field renormalization constants SZj; are general 2x2 matrices. Inserting fl3.65p and 
f l3.68p into the Lagrangian (13.641) . one obtains the counterterm Lagrangian 



= {fi 
- ( n n ) 



fi 

f2 



6M, 



Hence the renormalized self-energies for the sfermions are 

= ^f^^{p') + lp'(6Z} + 6Z^-U6Z}Mf + Mf6Zj 



For convenience we define 



5m| 

Ji 



Ur5MrU\) 

■' ■' J/ii 



5m 



/12 



Uf5MfU\ 

' ■> f/ 12 



(3.69) 



(3.70) 
(3.71) 

(3.72) 



The independent parameters in the sfermion sector are the soft-SUSY breaking parameters 
M? , M? , and At. Their counterterms are determined in the sfermion sector. Here we 
treat squarks and sleptons separately. 



3.6.1 Renormalization Constants for the Squarks 

There are five independent parameters: M? , M?^, M| , Au and A^ in every generation 
of squarks. In order to fix their counterterms, one can renormalize two up- and one of the 
down-type squarks via the on-shell renormalization scheme. Here we choose the lighter 
down-type squark di. The on-shell renormalization conditions can be written as follows, 

Retjjml) = 0, Retf-Jml)=0, i?eS)Jm|)=0, (3.73) 

where z = 1, 2 is the index of the squarks. / can be an up- or a down-quark. Inserting the 
renormalized self-energy (13.701) into the on-shell conditions above, and choosing 



5Z f = 5Z f , 

/12 /21 ' 



(3.74) 
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one gets the counterterms of the mass matrices and fields for the squarks. They can be 
written as 



ai 



5m 7 

J12 



SZ? 

ii2 



SZf 

J ii 



1 1 

m - — m - 

Jl /2 



-ReY!, (m| ) . 



(3.75a) 
(3.75b) 

(3.75c) 

(3.75d) 
(3.75e) 



The counterterms (5M? , SM? , 6Au are formulated from fl3.75al) and fl3.75cp {f = u), 
which can be expressed as 



1 



(3.76a) 
(3.76b) 



6Au = — 



Uui^Uu^2{5ml^-5ml^) + {Uu^,Uu^^ + Uu^^Uu2j5mi:,^^-5Cii,^ . (3.76c) 



The other countertems are derived from fl3.75bl) and fl3.75cl) {f = d). Their expressions are 



- U- U- f/- 

dl2 1^12 1^12 



(3.77a) 



Tjl TT- TT- 



J " Sm- - 2 

^ JJ2 "'"'U2 ^ fj2 



di 



di2 



di2 



di2 



SA. = — 



1 

rrid 



u. 



di "12 



+ 2 



dl2 ^dl2 
^dll^"12^«22 



(3.77b) 



£^12 



u. 



di2 



The counterterm 5mj can be expressed by the counterterms of the soft-breaking para- 



meters, 



(3.78) 



46 



CHAPTER 3. RENORMALIZATION OF THE MSSM 



3.6.2 Renormalization Constants for the Sleptons 

The independent parameters in each generation of the sleptons are M? , My , and Ai. Their 
counterterms are determined by renormahzing the two charged sleptons in the on-shell 
renormalization scheme. From the previous discussion, one gets the on-shell renormaliza- 
tion conditions for the slepton sector, 



Re±f (m? ) = , Re± 



0, 



i?eEf„(m?) = 0, Re%iml] 



0, 

: 0. 



(3.79) 



Similar to the squark renormalization, the counterterms and the field renormalization con- 
stants are obtained by solving the equations fl3.79p . 



6mf 



6mr 

112 

6Zr 

''ii 



6Zf 



''ii ^ H 

ReEj (m?) + i?eEr (m 



Hi H 

ReEf (mf)-ReEj (m? 



12^ 



-i?eS;-Ki). 



2 

m~ 

h 



(3.80a) 

(3.80b) 
(3.80c) 

(3.80d) 
(3.80e) 



From fl3.80al) and fl3.80bl) . one finds the expressions for the counterterms SM? ,SM? and 



SM? 



5Mf 



5Ai = — 
mi 



Uf 6mj + U? Smf - 2Uj Uj Sm, - 6Cf 

hi h li2 h '12 '22 '12 '11 

U? 6mj + Uf 6mf + 2Uj Uj 6mr - 6Cf 

ll2 h hi h '12 '22 '12 '22 

1 



f^^^iif^i2 (K - K) + ^UlriUl22 + Uh,U-,J5m-,^^ - SCj 



(3.81a) 
(3.81b) 

(3.81c) 



The counterterm for the sneutrino mass Mj.^ is not independent. It can be expressed by 
6Mf as in dSjlb]). 



3. 7. RENORMALIZATION OF THE HIGGS SECTOR 



47 



3.7 Renormalization of the Higgs sector 

In the CP-conserving MSSM the bihnear terms of the Lagrangian which describe the MSSM 
Higgs sector can be written as 



hi 
h2 



(3.82) 



where the Higgs multiplet {hi h2 ) can be ( A° ) , { h° H° ) oi { H"^ G^ ). 
Mj^ is the mass matrix of the Higgs bosons, Mj^ = UisM^oUj^, UaM^oU\ or Uf3M^±Ujj, 
respectively. Their definitions and expressions can be found in (12.651) and (I2.73p . 

In order to renormalize the Higgs sector we introduce the renormalization constants for 
their mass matrices and fields by the transformations 



M^o - 

hi 
h2 



1 + ^SZf, 



hi 
h2 



(3.83) 



where 6Zh is a general 2x2 matrix. Inserting these transformations into the Lagrangian 
fl3.82p . one can get the counterterm Lagrangian 



6C = ( h\ hi ) [-p' [6Zh + 6Zl 



M^6Z, + 6ZlM^ 



5M^ 



hi 
h2 



. (3.84) 



The renormalized self-energies for the Higgs bosons can be derived from the counterterm 
Lagrangian: 



S.,(/)+[/ {6Z + 6Z^)~-(m^6Z, 



. (3.85) 



The independent parameters in the Higgs sector are chosen to be the tadpoles T^o, Tfjo, 
the mass m^o, and the ratio of the vacuum expectation values tan/3. Their counterterms 
are fixed in the Higgs sector. 
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3.7.1 Counterterms for the Tadpoles 

The tadpoles T^o and Tj^o are equal to zero at tree-level. Their counterterms can be fixed 
by requiring that the renormalized tadpoles, which are defined as the unrenormalized ones 
plus the counterterms, are equal to zero at one-loop order as well, 

Refho = ReThO + 6ThO = , (3.86) 
Refffo = RcTho + 6Tho = . (3.87) 

Hence the counterterms for the Tadpoles can be expressed as 

5Th = -ReThO , 5Th = -ReTffo . (3.88) 



3.7.2 Counterterm for tan/3 



Since tan/5 is the ratio of the vacuum expectation values, tan/5 = ^, 
renormalization constants for the vacuum expectation values vi and V2, 

1 

Vl 



V2 



we introduce the 



1 + -5Zh, ) [vi + 5v{) 



1 + i^^Zh, 



(f 2 + 5V2 



(3.89) 



5Zhi and 5Zh2 are the field renormalization constants for the Higgs doublets Hi and H2, 
respectively. 



V2 

Vl 



(1 + \5ZH2) {V2 + SV2) one-loop V2 

{1 + 16Zh,){vi + 6vi) ~ 

Vl \ 2 V2 



l + \5Zu, + '-^ 
Vl 1 + \5Zh, + ^ 



-5Zh^ 



5vi 

Vl 



V2 
Vl 



1 + ^SZh^ - ^SZh, 



where we take ~ = 0. The transformation (13.901) can also be expressed as. 



S tan P = tan P 



SZh2 — SZhi 



The Lagrangian of the coupling of the neutral Higgs boson to Z boson is 

C^z = ifnz (p'^XiZf^ cos /? + p''x2^/. sin /?) , 
After renormalization, its counterterm Lagrangian can be written as 



^AOZcT 



-imzp'^A^Z^ sin P cos /? {5Zh2 — SZh^ 
imzP^C^Z^ 



,,^0 y 5Za^G^ + 5Zzz ^ ^^^^p,G'Z, . 



(3.90) 

(3.91) 
(3.92) 



(3.93) 
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Hence the renormalized self-energies are expressed as 

^AOz(/) = S^o^(p2) - imz sin p COS p (5 - 5Zh^] 
±coz{p') = ^cozii) + ^ ( + + 5mz ] . (3.94) 



According to [39] one can determine (5tan/3 by the requirement that the — Z mixing 
vanishes for an on-shell A^ boson, 

ReJ^Aozim^o) = ^ — — = — -Im ReT^Azim^) . (3.95) 

tan (j 2mz sm p cos p L J 



A convenient choice is the DR renormahzation of tan (3 [40] , which means that the coun- 
terterm only contains the UV-divergent parts, 

'^tan/jP^ 1 „ . 2^^ q«n 

tan jj Iniz sm p cos p 

Here the subscript "div' means that only the UV-divergent parts are considered. Since 
this choice has the advantage of providing the gauge invariant and process independent 
counterterms, it has been assessed to be the best choice of defining tan/3 [18]. 



3.7.3 Renormalization Constants for the Neutral CP-odd Higgs 
Bosons 

The renormalized self-energies for the CP-odd Higgs bosons can be expressed as 

Sao AO (p^) = Saoao(p^) + (p^ - ?T2,^o) 5^AOAo - 5m^o , (3.97a) 

Sgogo(p^) = Sgogo(p^) + {p^ - mlo) SZgogo - 5m^o , (3.97b) 

Saogo(p^) = Saogo(p^) + (^^AOGO + (^^GOAo) 

— -m\o6ZAOG0 — -niGodZG^A" — Sm\oGO ■ (3.97c) 



From fl2.73dp and fl2.73el) . one gets the explicit expressions for the counterterms Sm'^OGO 



and Sm^o, 



Smioco = (STho sin(a — (3) + 6Tho cos(a — /?)) , (3.98a) 

2mzSwCw 



SttIgo = (— 5Tf/o cos(a — /9) + 5T/jO sin(a — /5)) , (3.98b) 

2mzSwCw 

which are dependent on the counterterms ST^o and 6Tho. 

The counterterm for m^o is determined by renormalizing the neutral CP-odd Higgs 
boson A^ via the on-shell renormalization scheme, 

i?eSAOAo(m^o) = 0, (3.99) 
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which makes the renormahzed mass equals to the pole of the propagator. The diagonal en- 
tries of the field renormalization matrix are fixed such that the residues of the renormalized 
propagators are equal to 1, 



AO AO [''^ AO, 



, ReK 



GOGoyJ^GO) 



0. 



(3.100) 



The on-shell renormalization scheme also requires that the renormalized IPI two-point 
function for the CP-odd Higgs boson is diagonal for the external on-shell particles, which 
determine the non-diagonal entries of the field renormalization matrix, 



i?eSAOGo("^Ao) — 0' -ReS^OGo(m^o) = 0. 



(3.101) 



Inserting the expressions for the renormalized self-energies fl3.97p into the on-shell renor- 
malization conditions above, we obtain the expressions for the renormalization constants. 



5m\o 

6 Z AO AO 

5Zgqgo 

5 Z GO AO 



6Zaogo 



i?eS^o^o(m^o) , 
-ReE' 



AO AO AO) ; 

2 ( -ReT. AO GO im\o) + Sm\oGO 



m 



AO 



m 



GO 



2 f — _ReS^OGo(m^o) + 5m 



AO GO 



2 2 



(3.102a) 
(3.102b) 
(3.102c) 

(3.102d) 
(3.102e) 



3.7.4 Renormalization Constants for the Neutral CP-even Higgs 
Bosons 

The renormalized self-energies for the CP-even Higgs bosons can be expressed as 



^HOHO {p 



^hORO 



ShOho(p2) + [p^ - mlo) SZhOho - 5mlo , 
^HOHo{p'^) + [p^ - m\jo) 6Zhoho - Smjjo , 

^hOHo{p^) + \p^ {SZhOHO + SZhoho) - 



[m\o6ZhOHO + m]joSZHOho) — Sm 



2 

hORO 



(3.103) 
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From fl2.73ap . fl2.73bl) and fl2.73cl) one derives the expressions for the counterterms Smf^o, 
6mjjo and 5m\o^o, 



5m\o = Sm\o cos^(q; — j3) + 5m\ sin^(Q; + /?) + 

-5Tho cos(q; — (3) sin^(a — (3) + 



2mzSwCw 

5Tho sin(a — + cos^(a — 13)) + 



2mzSwCw 

5 tan (3 cos^ (3 {m\o sin 2 (a — (3) + m| sin 2(a + (3)) , (3.104a) 
5m\o = 6m\o sin^(a — /?) + 5m| cos^(q; + /?) — 

5Tho cosia — P) (l + sin^fa — P)) — 

2mzSwCw 

-(5T/jO sin(a — /?) cos^(q; — /?) — 



2mzSwCw 

S tan P cos^ /? {m\o sin 2 (a — p) + m| sin 2(a + p)) , (3.104b) 
Smloffo = - ((5m^o sin2(Q; — /5) — 5m|sin2(a + /5)) + 

{6Tho sin^a - P) - 6Tho cos^(a - p)) - 



2mzSwCw 

(5tan/3cos^/3(m^oCos2(a-/5) + m|cos2(a + /3)) , (3.104c) 



where all the counterterms appearing on right-hand side are already determined. Hence 
in the neutral CP-even Higgs sector, only the field renormalization constants are to be 
fixed. Similarly to case of the neutral CP-odd Higgs boson, the on-shell renormalization 
conditions for the CP-even Higgs sector can be formulated as 



Ret'^Oho{mlo) = , Ret'HOHoirriHo) = , 

i?eEfto/i-o(m^o) = 0, i?eS/iO/^o(m^o) = 0. (3.105) 



Applying the expressions for the renormalized self-energies (13.1031) to the on-shell condi- 
tions above, one obtains the field renormalization constants. 



-ReT,hom(mln) + 5mln^a 
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3.7.5 Renormalization Constants for the Charged Higgs Bosons 

The renormalized self-energies for the charged Higgs bosons can be expressed as 

^H-H+{p^) = T.H-H+ip'^) + [p^ - m\j±) 5Zh-h+ - 5mjj± , (3.107a) 
^g-gAp') = ^G~GAp') + {p'-^l±)SZG-G+-6ml^, (3.107b) 

^H-G+iP^) = ^H~G+{.P^) + ]^P^ {5Zh-g+ + 5Zh-g+) 

- ^m%±6ZH-G+ - ]^m%+5ZG+H- - 5"^h-g+ • (3.107c) 



From (]2.73f|) one gets the exphcit expressions for the counterterms Sm'jj±, (5m^_^+ and 

5m'\j± = 5m\o + Sm^ , (3.108a) 

5m?r-f^+ = {6Tfjo sm{a — (3) + 6ThO cos{a — (3)) 

2mzSwCw 

- 6tan(3mjj±cos^ p , (3.108b) 

Smf.± = (—6Thocos(a — p) + SThosm(a — P)). (3.108c) 

2mzSwCw 

The field renormalization constants are determined by the on-shell conditions, 

i?eE'^-^+(m^±) = =^ Zh-h+ = ~ReT,'jj-jj+{mjj±) , (3.109a) 
Re%-GAml±) = =^ Zg-g+ = -ReJ:'G-GAml±) , (3.109b) 

ReJ:H~G+{ml±) = =^ 5Zg+h-=2 " \ ^^-^ , 3.109c 

ReT.H-G+irnl±) = Q =^ 5Zh-g+='2^ ^/ ^ , ^^-^ . 3.109d 

One can also fix the field renormalization constants for the Higgs fields via the DR scheme, 
where the counterterms contain only the UV-divergent parts, 

5Z^ = . (3.110) 



Chapter 4 



Calculations for the Next-to-lightest 



The Minimal Supersymmetric Standard Model (MSSM) is the supersymmetric extension 
of the Standard Model (SM) with the minimal particle content, the details have been 
presented in Chapter [2l In the MSSM with conserved i?-parity, the lightest supersymmetric 
particle (LSP), which in many scenarios is the lightest neutralino Xi, appears at the end of 
the decay chain of each supersymmetric particle. The LSP escapes the detector, giving the 
characteristic SUSY signature of missing energy. While this helps to suppress backgrounds 
from SM processes, it also makes the measurement of supersymmetric particle masses more 
difficult. 

At the LHC, the total SUSY production- cross section is expected to be dominated by 
the production of gluinos and squarks, which decay into lighter charginos or neutralinos. 
Of particular interest are decay chains leading to the next-to-lightest neutralino X2- X2 
turn can always decay into the LSP and two fermions //, at least for light SM fermions 
/. The leptonic final states are of particular interest, since they can be identified relatively 
easily even at the LHC. Depending on neutralino, slepton and Higgs boson masses, the 
possible leptonic decays of X2 are three-body decays X2 ~^ Xi^~i~^y cascade two-body decays 
X2 ^^^^ Xi^~^~^ and/or X2 ~^ ~^ Xi^~^~^j where 0° stands for one of the three 

neutral Higgs bosons or the neutral Goldstone boson of the MSSM. The dilepton invariant- 
mass distribution of these decays has a specific shape with a sharp edge near the endpoint 
which only depends on the kinematics. This distribution can be measured accurately at 
the LHC. In particular, its endpoint is used in several analyses that aim to reconstruct 
(differences of) supersymmetric particle masses [H E]. Under favorable circumstances it 
has been shown that the endpoint can be measured to an accuracy of 0.1% at the LHC 
In order to match this accuracy in the theoretical prediction, at least one-loop corrections 
to X2 decays have to be included. 

Turning to the planned e^e~ linear collider ILC, X1X2 production is often the first 
process that is kinematically accessible [19] (other than x^ pair production, which leads to 
an invisible final state). The detailed analysis of X2 decays can then yield information about 



Neutralino 




Decay 



53 



54 



CHAPTER 4. CALCULATIONS FOR xl DECAY 




Figure 4.1: The Feynman diagram for the squark decay Qs qX2 ~^ IXi^ cit the LHC, 
s = 1, 2 labels the squark mass eigenstates. 



heavier supersymmetric particles. Under favorable circumstances, (9(10^) X2 Xi^~^^~ 
decays may be observed at the ILC, again making the inclusion of quantum corrections 
mandatory to match the experimental precision. 

In this chapter, we calculate leptonic X2 decays at one-loop level. Cases where X2 has 
two-body decays X2 ~^ ~^ Xi^^^^ (^1 stands for the lighter one of the two charged 
sleptons) are treated both completely and in a single-pole approximation. In the complete 
calculation one has to employ complex slepton masses in the relevant propagators and 
one-loop integrals. The single-pole approximation in this case is performed in the way that 
the X2 decays are treated as the production and decay of the sleptons Zi. We compare the 
results from the complete and approximate calculations and find a good agreement. We also 
analyze a scenario where X2 only has three-body decays. The virtual photonic contributions 
are infrared (IR) divergent, hence the contributions of the real photon bremsstrahlung must 
be added to the one-loop corrrections in order to cancel these IR divergences. In addition 
to calculating the integrated partial widths, we study the differential decay width of X2 
a function of the dilepton invariant mass. In order to obtain the total decay width of X2 
and hence the branching ratios of the leptonic decays, the invisible decays X2 ^ Xi'^i^i and 
the hadronic decays X2 ~^ Xil'l are also calculated. 

This chapter is organized as follows. Section 14.11 gives the tree-level calculations for 
X2 leptonic decays. In Section 14.21 we discuss how to calculate these decays completely at 
one-loop level, where the virtual corrections and the real photon bremsstrahlung are con- 
sidered in detail. When the lighter sleptons Zi can be on shell, these decays are calculated 
approximately in Section 14. 3[ The total decay width of X2 ^^^d the branching ratios of the 
leptonic decays are discussed in Section 14.41 where the invisible decays and the hadronic 
decays are calculated. The numerical results and discussions are given in Section where 
the SPSla parameter set [71 [8] is presented in detail. 



4.1. TREE-LEVEL CALCULATIONS FOR ^ X^L'L^ 

4.1 Tree- level Calculations for X2 ~^ Xi^^^ 
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The Bom Feynman diagrams for X2i^ 



Xi{k2)l {ks)l^{ki){l = e,/i, r) are displayed 



in Figure I4.2[ The Mandelstam variables are defined as 

T = (ki- k^f , T12 = {ki - k2f , 5*23 = {k2 + ^3)^ 
The decay width of this process can be written as (see Appendix IC.2 



'1^3 - 



(4.1) 



(4.2) 



where M*^°) is the matrix element of the Born diagrams, it is squared and averaged over 
the spin of the external particles. The expressions for the phase-space element (i$i^3 can 
be found in Appendix I C. 2 [ 

If the two-body decays X2 ~^ ^1^^ ~^ Xi^^^^ are kinematically allowed, i.e. the sleptons 
li can be on shell at some points in the phase space, a finite width of li is necessary. It 
arises from the imaginary part of the slepton self-energy. A finite width is introduced via 
Dyson summation 



+ 



h h 



/j2 _ ^2 



+ ■ 



F - m? + S(fc2 

'1 



(4.3) 



where S(/c^) is the renormalized h self-energy. Following Ref.[53j, a gauge invariant matrix 
element is obtained by a Laurent expansion around the complex pole: 



1 



P - m? + S(A;2 



— 



ReS(P) 
k'^ — rn^ 

h 



(4.4) 



were mi denotes the position of the complex pole in 04.31) . It is obtained as the solution of 



(4.5) 




Figure 4.2: The Born Feynman diagrams for X2 — ^ Xi^~^'''(^ = e, /x, r), s = 1, 2 labels the 
slepton mass eigenstates, denotes the MSSM neutral Higgs boson /i", A^ and neutral 
Goldstone boson . Since the Yukawa coupling (f)^l~l~^ is proportional to the lepton mass, 
the Higgs intermediate states are neglected when / = e and /i. 
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For the tree-level amplitude the complex pole is calculated at one-loop level. Its 
explicit expression is 

nip = m? — imii_^T^^^ , (4.6) 

where we have employed on-shell renormalization scheme as in Section 13.6.21 and T~^^ is 

the tree-level decay width of ^i, mi_Tj^^ is the imaginary part of the slepton self-energy 
S(m? ). Since the second term in the parentheses in fl4.4p is at one-loop level, we do not 

^ 1 ^ 

need it in the tree-level calculations. Therefore, the gauge-invariant tree-level amplitude 
for the decays Xo, ~^ Xi^~^~^ can be written as 

= ^^ifiil__Azi!il_^ + , (4.7) 

— rrij + imr Tr 

where V^-±^^ and ^-±^o;± represent the X2^iP and ifxil^ vertices, respectively, B{k'^) 
denotes the non-resonant part of the matrix element, i.e. the matrix element of the diagram 
(a) and (b) for s = 2 and diagram (c) and (d) in Figure IT2l 

Using the residue theorem one can easily obtain the relation for an analytic function 

f . -^Z^AO- itrf«m,. (4.8) 



_ ^2 _|_ jj^K )^ 
h h '1 



mr T- 
'1 h 



This means that the function f{k'^) in (14. 8 p will be dominated by the regions of k"^ close 
to m? if vf^ ^ nij . The non-resonant part of the matrix element is much smaller 
than the resonant one (diagram (a) and (b) for s = 1 in Figure 14.21) . hence it can be 
neglected approximately. Applying (14. 7[ 14.81) to (14.21) . neglecting the contributions from 
the non-resonant diagrams, the decay width of X2 ~^ Xi^~^~^ can be written approximately 
as 



''X2 

1 1 



(2vr)2 2m,-^rf y ^' 

= rW(x°-/fni?r(°)(/^^X?/^), (4.9) 
where the branching ratio of the decay if — > Xi^^ is defined by 

Br^'\!f^Xr) = ^^-^^i^^. (4.10) 

Hence, when the lighter sleptons li can be on shell, we can compute the relevant partial 
widths in the single-pole approximation, where the decays X2 ^ Xi^~^~^ ^-re treated as the 
production and decay of the sleptons h. 
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4.2 Complete One- loop Calculation for X2 ~^ Xi^^^ 
4.2.1 Virtual Corrections 

In general the virtual one-loop corrections to three-body decays can be classified as vertex 
contributions, self-energy contributions and box contributions. The first two classes are 
UV finite after adding the contributions from the counterterms that originate from the 
renormalization of the MSSM, as discussed in Chapter [31 The box diagrams are by them- 
selves UV finite. Different types of diagrams and their counterterm diagrams are shown 
in Figure (14.31) . The MSSM Feynman rules, as well as the resulting counterterms, are im- 




Figure 4.3: Classification of the one-loop virtual diagrams 



plemented in the FeynArts package of computer program ^Slj , which allows an automated 
generation of the Feynman diagrams. The matrix element and the one-loop integrals are 
calculated with the help of the packages FormCalc and LoopTools [52], respectively. The 
squared matrix element at one-loop level can be written as 

|mW|' = |m(o) + M,irt|' ~ |m(°)|' + 2Re[M(°)Mtrt] , (4.11) 

where M^-„t is matrix element of the virtual one- loop diagrams for the process X2 ~^ Xi^~^^- 
The generic virtual one-loop diagrams are shown in Appendix Applying (14.111) to the 
expressions for the width (14.21) one obtains 

r(i) = r(°) + rvirt. (4.12) 

Similarly to the tree-level case, diagrams with a slepton li propagator have singularities 
when li can be on shell. The single-pole approximation can also be used at the one-loop 
level, which will be discussed in Section 14.31 Here we focus on the complete calculation. 
Following the strategy in Section 14.11 one can obtain a gauge invariant matrix element 
at one-loop level. In order to obtain 0{a) accuracy near the h resonance, one needs to 
calculate the complex pole to two- loop level [53] . 
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where we have apphed the on-shell renormahzation scheme at two-loop level, and r[^'' 

denotes the one- loop-level width of li. Then the gauge invariant matrix element at one- 
loop level can be written as 



h '1 h 



(4.14) 



where C{k'^) denotes the non-resonant part of the matrix element, the residue 74(/c^) can 
be expressed as 



A(k^ 



ReS(A;2 



V. 



.2\T>(1) 



(1) 



h 



(4.15) 



where V^lj±i^^ and ^-±^o;± represent the renormalized X2^fl^ ^^'i ifxi^^ vertices at one-loop 
level, respectively. 




(d) 



Figure 4.4: Virtual photonic corrections in the complete calculation, s = 1,2 labels the 
slepton mass eigenstates, (j)^ denotes the neutral Higgs boson h^, H^, and the neutral 
Goldstone boson 
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Moreover, the one-loop integrals in the diagrams shown in the first two lines of Figure|43] 
also give singularities when the sleptons li are on shell. One should therefore use complex 
slepton masses, 

^ ^ (4.16) 



- ml k^-m^ + imr 

in the one-loop integrals from these diagrams. The one-loop integrals with complex masses 
can be calculated automatically by LoopTools. The analytical expressions for scalar three- 
point and four-point functions with real arguments can be found in Refs. pTilMt l55| l56] . We 
generalized the scalar four-point function to allow for complex arguments. The scalar three- 
point functions from the diagrams (al), (a2), (bl) and (b2) in Figure are calculated 
analytically. The analytical results are presented in Appendix [Bl 

4.2.2 Analytical Results for Virtual Photonic Corrections (/ = e,/x) 

The virtual photonic diagrams are shown in Figure where the diagrams (a3), (b3), 
(c) and (d) have the property that the virtual photons are attached to external on-shell 
charged particles. This results in IR divergences, which we regularized by introducing 
a fictitious photon mass A. The IR divergences cancel after we add contributions from 
real photon bremsstrahlung, which will be discussed in Section fl4.2.3p . The masses of the 
light leptons, i.e. mi (/ = e, /i), are neglected except when they appear in the one- loop 
integrals. The virtual photonic corrections are calculated analytically, where the terms 
which contain the soft singularity In A, the mass singularity Inm; and the large logarithm 
ln(mj- /r^-) (only exits when li can be on shell), are treated separately as the singular part 
of the virtual photonic corrections. 

The photonic corrections from the diagrams in the first two lines of Figure 14.41 (with 
their counterterms) can be written as 

= ^Q?(wEl^i°^l' + ^-,.$:|Mrr)rf$i^3, (4.17) 

where Mi°^ and M^^^ are the matrix element of the Born diagram (a) and (b) in Figure 
14. 2^ respectively, Qi denotes the lepton charge, Qf = 1. The factor (5virt,a can be expressed 
as 



-'virt,a rn 9 

I — m^o 

,2 



(So(m2o,m2,m2-) -5o(r,0,m2-)) + 

/ \ 
(Bo{mlo,mlml) - Bo{T,0,ml)) + 

X? 

2 (^BoimlmlO) - ^5o(T,0,mJ)^ - (mlo - mfJC^ - 

(m^o - mfJC', - TuiT - mfjDo + 2 SZ^l , (4.18) 



T 



60 



CHAPTER 4. CALCULATIONS FOR xl DECAY 



where Cq, Cq and Do are scalar three- and four-point functions, 5Zj is the photonic part 
of the lepton field renormalization constant SZ^f^ or SZll^ which have been expressed in 
(]3.25bl) and fl3.25cp . respectively. The general definitions of the scalar one-loop integrals 
can be found in Appendix [Bj The arguments of Cq, Cq and Dq in fl4.18p are expressed as 
follows, 



Co = Co(mf , m^o, T, 0, , mfj , 
Co = Co{mf, m^o, T, 0, mf, mfJ , 

Do = Do{mf,m\,m\,rnf,T,Ti2,0,m'f,mf^,m'f) . (4.19) 



As discussed before, one should use complex masses m? 



m? +11717 rp'' in the scalar three- 
is h 

and four-point functions shown in fl4.19p when the slepton 1^ can be on shell. The scalar 
two-, three- and four-point functions as well as the lepton field renormalization constant 
in (I4.18P are calculated in Appendix [B], where the three- and four-point functions with 
complex arguments are also presented. One finds that the expressions in fl4.18p are UV 
divergent. This indicates that the virtual photonic corrections are UV divergent. Adding 
these UV-divergent part to the non-photonic virtual corrections, one obtains UV-finite 
results. 

The terms in ( 14. ISp which contain the soft singularity In A, the mass singularity Inm; 
and the large logarithm ln(m;-/r;-) are treated separately as the singular part of ^virt.a, 
which can be expressed as 




In 



mt 



+ 



(4.20) 



where the terms proportional to lnr(m^ — T — iT^-^m,-) /mf] only exist when the slepton 

Is can be on shell. The factor 6vht,b in (I4.17P is calculated in the same way as 5virt,a- Its 
singular part can be obtained by replacing T with S23 in (14.201) . 



"virt,fe 



virt,a 



(4.21) 



The photonic virtual diagram (c) in Figure 14.41 is calculated and its contribution can 
be written as 



dT 



virt 



1^3 



(4.22) 
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where Mc^^ is the matrix element of the Born diagram (c) in Figure 14.21 The factor 5virt,c 
reads 

5virt,c = -Qri2 + 25oK,0,m2)-^5o(Ti2,m2,mn + 2 5Zj, (4.23) 

where Cq is the scalar three-point function Cq = Co{m'f,m'f,Ti2,0,m'f,m'f). The scalar 
two-point and three-point functions as well as the lepton field renormalization constant in 
f l4.23p are calculated in Appendix [Bl One finds that the factor 5virt,c is UV finite. The 
singular part of 5virt,c; which contains the singularities In A and Inm;, are 

Combining the expressions for the photonic virtual corrections (14.171) . (I4.20p . (14.211) . 
(I4.22P and (I4.24p . the singular part of virtual photonic corrections {I = e, /i) can be written 
as 



dr 



a 



virt 



smg 



TT 



(5: 



sing _ rsing ■ 
virt, a virt,c. 



virt,b 



+6"- 



virt,c 



(4.25) 



Here we have concentrated on the case where the lighter slepton li can be on shell. The 
factor S^^f^ — ^virf c proportional to In \{mf — T — iTf^mr ) / mf 1 . After integration over 
the Mandelstam variable T, one obtains 




„ \dT = iV^pmr In ^ I + ■ ■ • , (4.26) 

where we have set T = mf because in this region the integrated result is dominant. The 
ellipses in (I4.26P represent terms which have nothing with Tf\ Since 

^ 1 



lim iT^rrif In I ^ I = , (4.27) 

(1) «i m,- ' ^ ' 



-p(l) 

-li r 



the large logarithm In ^r|-^^/m^- j disappears in the integrated result when v'^p 0. Sim- 
ilarly there is a large logarithm In {v^ /rriiA in the factor ^Xfb ~ Cirf,c (14.251) . which 



also disappears in the integrated result when vf'' goes to 0. 

^ 1 

Since we keep the r mass everywhere, the virtual photonic corrections with r^r+ final 
states are calculated numerically. 



62 



CHAPTER 4. CALCULATIONS FOR xl DECAY 



4.2.3 Real Photon Bremsstrahlung 

In order to cancel the IR divergences in the virtual corrections, we have to add contributions 
from real photon bremsstrahlung to the one-loop corrections. The diagrams for the process 
X2(^i) — ^ Xi{k2)l~ {k3)l~^ {k4^)j{q) (/ = e,fi,T) are displayed in Figure 1^31 Generally the 
decay width of the real photon bremsstrahlung can be written as 

1 /" ^ 2 

Tbrems = TTTTsT^ / / .M^brems (i$i^4 , (4.28) 

(27r)«2m^o J 

where Mbrems denotes the matrix element of the diagrams in Figure 14.51 The definition for 
the phase-space element (i$i^4 can be found in Appendix IC.2[ One must use the complex 
slepton masses in the propagators as (14.161) when the lighter sleptons li can be on shell. 

From Figure we know that the real photon can be emitted from the charged leptons 
and sleptons. In the case of photon emission from a charged lepton, the amplitude for the 
real photon emission has the property 




This will give rise to IR divergences which are regularized by a photon mass. The "soft 
photon bremsstrahlung" contributions are defined via the condition that the photon en- 
ergy < AE, where the cutoff parameter AE should be small compared to the rel- 
evant physical energy scale. The complementary contributions are called "hard photon 
bremsstrahlung", which are defined as real emission contributions with E^ > AE. The 
contributions of soft photon bremsstrahlung are IR divergent, which are sufficient to can- 
cel the IR divergences in the virtual corrections, while the contributions of hard photon 
bremsstrahlung are IR finite. The contributions of the real emission can be expressed as 

Tbrems Tsoft 

(AE) + rhard(A£;). (4.30) 

The dependence on the largely arbitrary parameter AE cancels after summing soft and 
hard contributions, provided it is sufficiently small. The contributions of the photon radi- 
ation from charged sleptons are always finite because the sleptons are internal particles. 
If we neglect the light lepton masses, the factor ^ in Km can be written as 

(4.31) 



2ki ■ q 2kioE^{l — cos 6i^) ' 

where 6i.y is the angle between momentum of the photon and momentum of the emitting 
particle. If cos^/^ ~ 1, i.e. the photon and its emitter are collinear, the contributions of 
the hard photon bremsstrahlung also contain a divergence. It is regularized by the masses 
of the leptons in the final state. However, since the lepton masses, i.e. rUe and m^, are very 
small, it is very difficult to get stable numerical results from a direct numerical evaluation 
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I- l- 




Figure 4.5: The Feynman diagrams for the real photon bremsstrahlung '}^{ki) — > 
Xi{k2)i~ {k3)l~^ {k4)j{q) , s = 1,2 labels the slepton mass eigenstates, 0" denotes the neutral 
Higgs boson H^, A^ and the neutral Goldston boson G^. 

of hard photon bremsstrahlung even we keep the light lepton masses everywhere. This can 
be overcome by dividing hard photon bremsstrahlung into a collinear part, where the angle 
between the photon and the radiating particle is smaller than a very small angle A^, and 
the complementary non-collinear part, 

rhard(A^) = rcoll(AE,A^^)+r„on-coll(AE,A^^). (4.32) 

The angular cutoff A9 should be so small that we can assume that the collinear-photon 
emission does not change the direction of the three-momentum of the emitting lepton. 

So far we have divided the four-particle phase space into a soft, a collinear and a 
non-collinear region. This is called phase-space-slicing method. The soft and collinear 
contributions can be calculated analytically, while the non-collinear contributions are cal- 
culated numerically using the multi-channel-approach in the Monte Carlo program [571 EH] • 
The details about this technique are presented in Appendix O 
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Soft Photon Bremsstrahlung 

Since the energy of the emitted soft photon is by definition very small, this emission 
essentially does not change the momenta of the other final state particles. In the soft 
region, the squared amplitude |Msoft| can be written as the Born squared amplitude 
|M'^°)|^ multiplied by a soft factor. 



soft 



2„2 



-Qte 



^3 



3-q 



C4 ■ q 



(4.33) 



where we have neglected all the terms proportional to the photon momentum q in the 
numerator. The four-particle phase space can also be factorized into a three-particle phase 
space and a soft part. Therefore the contributions of the soft photon bremsstrahlung can 
be written as 



dV 



soft 



a 



27r2 
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2^ 



d'\q\ 



^4 



Q'(/33 + /44 - 2J34)|M(°f ^$1. 



1^3 



(4.34) 



where 



E. 



hj 



d'\q\ 



\q\<AE 

2E^ {h-q){kyq) ' 



(4.35) 



7 - V lyi +^'- 

The general analytical expressions for the integrals lij can be found in Refs. [11], [51] . 
This has been implemented in the package of FormCalc. The light lepton mass, i.e. rrti {I = 
e, n) is kept only when it acts as a regulator for the mass singularity, hence the analytical 
expressions for the soft contributions can be written as 



dT 
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(4.36) 



where fcjo denotes the energy of the charged lepton whose four- momentum is defined as fcj, 
the dilogarithm Li2{x) is defined in Appendix iBl 
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Collinear Photon Bremsstrahlung 

From the discussions above we define the coUinear region: > AE and l—6g < \ cos^^^l < 
1, where 6g = 1 — cos A6'. This means that the colhnear part describes real photon radiation 
outside the soft photon region and cohinear to the emitter. We consider a final-state 
radiation, 

l{ki) liq) 2 = 3,4 (see Figure USD- (4.37) 
The squared matrix element of the collinear photon bremsstrahlung can be written as 

where Zi = P{zi) = jz^, ki = ^. The phase space in the collinear region can be 

separated into a three-particle phase space multiplied by a collinear factor, 

= — — ; — z o^'iki — ko — kr^ — kA) 

2k2o 2k3o 2A:4o 2E^ ^ ^ ^' 

= d^._^.^c^,z^^. (4.39) 



This is corresponding to the phase space of the radiation ^(^3) 7(g) + l{k^). Therefore 
the differential collinear contributions can be expressed as 

4 

t^reoii = ^Q] Yl dr^'\k,), (4.40) 

j=3 



where 



G^ = - T—[Piz^)'Tr^ Z (4.41) 

IT J ki-q \ ki-qj 2E^ 



It can be written as 



0/ M I 4fcjo 2^ 2zi 
P[Zi) In — 5-77 ~ ~^ 

ml 2 ' 1 - z,. 



dzi (4.42) 



after integrating out the solid angle of the momentum q analytically. 

If we treat a charged lepton and a collinear photon inclusively, i.e. the momentum 
of collinear photon is added to that of emitting lepton, the variable Zi in fl4.42p can be 
integrated out analytically. Hence the differential contributions of the collinear emissions 
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are written as the differential tree-level decay width multiplied by a universal function 
[591 EOj, 



ZTT ^ — ' 



3 , /AE\ 
- + 2 In ^ 

2 V ; 
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(4.43) 



This approach is for collinear-safe observables [60]. If one adds the soft and collinear 
contributions to the virtual corrections, all singularities (Inm^ and In A) cancel. This is 
in accordance with the Kinoshita-Lee-Nauenberg theorem [HI]. At the LHC the electron 
energy is determined calorimetrically. In this case a collinear photon would hit the same 
cell of the calorimeter as the electron, so the two energies cannot be disentangled. Hence 
the electron observables are defined as collinear-safe observables in our calculation. 



We also consider the non-coUinear-safe observables [60], i.e. the lepton and its collinear 
photon are not treated inclusively. Since the phase space depends on the variables Zi, the 
integration over Zi cannot be performed analytically. The differential contributions of the 
collinear photon bremsstrahlung are written as 



dr 
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(4.44) 



In this case the mass singularity In m; cannot be canceled in the differential width and hence 
become visible. At the LHC, muons pass through the calorimeter, where the photons are 
detected, and measured forther outside in the muon detector. Hence the muon observables 
are treated as non-collinear-safe observables in our calculation. 



4.2.4 QED Corrections 



The virtual photonic corrections by themselves are UV divergent, hence one cannot mean- 
ingfully separate the QED corrections from the one-loop contributions by simply selecting 
diagrams which contain a photon. Since the light lepton (e and //) final states and the r fi- 
nal states are treated differently, the corresponding QED corrections are defined differently 
as following. 

In the case of the light lepton final states, the photonic virtual corrections and the soft 
photon bremsstrahlung are calculated analytically. We can pick out potentially large QED 
terms from the the sum of virtual and soft photon bremsstrahlung corrections, dT^^^t+dTsoft, 



dT^ij-t + dVsoit 



(4.45) 



4.2. COMPLETE ONE-LOOP CALCULATION FOR X2 X{L'L 



67 



Here dV contains all the potentially large QED terms 
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(4.46) 



where we have assumed that the lighter slepton li can be on shell. The matrix element 



Mi°) and Ml""' are obtained from the Born diagrams (a) and (b) in Figure respectively. 



r(0) 



As discussed in Section 14.2.31 the terms proportional to In [v'^^^ /mi^^ in fl4.46p disappear 

in the integrated width, drremainder defined in f l4.45p is IR and UV finite and free of such 
large QED logarithms in (14.461) . The "QED contributions" can then be defined as 



dV 



QED 



dV + (iFhard • 



(4.47) 



Note that the QED corrections defined in this way do not depend on the cutoff parame- 
ters AE and AO. Moreover, terms proportional to logm^ in (14.471) cancel in the integrated 
width and in the differential width for the collinear-safe observables. Using the definitions 
(I4.45P and (I4.47p . the complete one-loop contributions can be written as 

= C?rW + rfr™nder + f^rQED • (4.48) 

When T~T^ are the final states of X2 decay, the r mass m,- is kept everywhere. This 
mass is so large that a stable numerical result can be obtained from the hard photon 
bremsstrahlung, hence we do not need to divide the hard photon bremsstrahlung contri- 
bution into coUinear and non-collinear parts. We follow a slightly different procedure to 
define the "QED part" of the correction. The virtual corrections contain photonic and 
non-photonic contributions. 



"-^ virt ^ "-^ virt 5 



(4.49) 



both of which are UV divergent, while the sum is finite. The photonic virtual corrections 
can be split into an UV-finite part dV and an UV-divergent part (iFuv-div' 



(4.50) 
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Here drjjy_^^_^ contains the terms that would be subtracted in the dimensional reduction 
regularization of rfF^j^- After this rearrangement, the virtual corrections can be written as 

= dT + rfFremaindcr 5 (4-51) 

where (iF remainder = dV^y^^^^ + dT^°^~'^ as well as dT are UV finite. The "QED corrections" 
are finally defined as 

dTQED = df + (ilbrems , (4.52) 

where (iF brems stands for the contribution from all diagrams with real photon emission. By 
construction, c^Fqed is UV and IR finite. 



4.3 Approximate One- loop Calculation for X2 ~^ Xi^ 

If X2 ~^ ^1^^ two-body decays are allowed and X2 does not have other two-body decay 
modes, at one-loop level, just like at tree-level, the decays X2 ~^ Xi^^^^ = r) can 
be approximately treated as production and decays of the slepton if, 

T^'Hx'^^xm ^ F«(x°-/^n5r«(/f-X?^^), (4.53) 

with 

5r«(/f^xr) = LJfL_^. (4.54) 

The virtual contributions of the production and decay of the slepton if, which now only 
contain vertex type corrections, are again calculated with the help of the programs Fey- 
nArts, FormCalc and LoopTools. The virtual photonic diagrams are shown in Figure IT6l 
which are IR and UV divergent. The scalar three-point functions corresponding to these 
diagrams are expressed analytically in Appendix [B] in the case of Z = e, /x. An IR fi- 
nite result is obtained by adding the contributions from the real photon bremsstrahlung 
(see Figure 14. 7p . where both the diagrams (a2) and (b2) must be included in order to 
preserve the gauge invariance. The contributions of the real photon bremsstrahlung are 
again separated into an IR-divergent soft part and an IR- finite hard part. The division of 
the hard photon bremsstrahlung contribution into a collinear part, which can be calculated 
analytically, and a non-collinear part, which is calculated numerically, proceeds along the 
lines described in Section 14.2.31 As discussed in Section 14.2.41 the UV-divergent photonic 
contributions cannot be treated separately as "QED corrections". We define the "QED 
corrections" in the same way as the complete calculation in Section 14.2.41 One finally 
arrives at a total one-loop contribution which is independent of the cutoff parameters. 
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Figure 4.6: Virtual photonic corrections in the approximate calculation. 




(al) (a2) (bl) (b2) 



Figure 4.7: Diagrams for the real photon bremsstrahlung in the approximate calculation. 

4.4 Total Decay Width of X2 ^^id the Branching Ratios 
of the Decays X2 ~^ Xi^^^~ 

As discussed in the beginning of this chapter, the next-to-lightest neutralino X2 decay 
into the LSP Xi and two fermions //. The leptonic final states are important because 
they can be identified at the LHC Moreover, the endpoint of the dilepton invariant-mass 
distribution is used to determine the mass relations of supersymmetric particles. The 
invisible X2 decay modes, i.e. X2 ~^ Xi^i^h do not effect the dilepton invariant-mass 
distribution. But they contribute to the total width of Xt Since it is very difficult to 
identify quarks at the LHC, the hadronic decays X2 ~^ XiQQ ai'G l^ss interesting than 
leptonic decays. In order to obtain the total decay width of X25 these hadronic decays 
must be calculated . The total decay width of X2 can be written as 

= E - l^l^xl) + ml - mx'i)] + E - qqx\) . (4.55) 

l=e,ii,T q=u,d,c,s,b 

Here we assume that the decay X2 ^ Xi^t is not kinematically allowed. The branching 
ratios of the leptonic decays X2 ~^ Xi^^^^ are defined as 

Brixl - xrn = ~Z . (4.56) 
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4.4.1 The Invisible Decays X2 xi^i^i 

The invisible decays X2 ~^ Xi'^i'^i ^-i'^ calculated at tree and one-loop level. The Born 
Feynman diagrams are shown in Figure 14.81 Here we focus on the case where the decays 




Figure 4.8: The Born Feynman diagrams for X2 — ^ Xi'^i'^i- 

X2 ~^ Xi^i'^i ell's pure three-body decays. Similarly to the calculations for the leptonic 
three-body decays X2 — > Xil^l^, these decays are also calculated with the help of FeynArts, 
FormCalc and LoopTools. The one-loop corrections are again classified as vertex, self, and 
box contributions. Since none of the external particles carries electric charge, there are 
no corrections involving real or virtual photons, and hence no IR divergences. Therefore, 
there are also no QED corrections in these decays. This makes the calculations for the 
invisible decays much simpler than for the leptonic decays. 

4.4.2 The Hadronic Decays X2 ~^ XiQQ {q ^) 

The hadronic decays of X2 calculated in order to obtain the total width of Xt The 
Born Feynman diagrams for the decays X2 ~^ XiQQ {q 7^ t) shown in Figure 14.91 where 
the masses of the light quarks, i.e. u, d and s, are neglected. Here we only consider the case 




Figure 4.9: The Born Feynman diagrams for X2 — ^ xHlQ {1 ^ s = 1,2 labels the squark 
mass eigenstates, 0° denotes the neutral Higgs boson /i^, A^ and the neutral Goldstone 
boson G^. 

where the hadronic decays X2 ~^ XilQ. ^"^^ pure three-body decays. Since the SUSY-QCD 
corrections are not considered in our calculations, the decays X2 ~^ XiQQ can be treated in 
the same way as X2 ~^ Xi^~^~^- 

The one-loop virtual corrections for the decays X2 ~^ XiQQ ^-^'s also classified into three 
types: vertex contributions, self-energy contributions and box contributions, where the 
vertex and self-energy contributions have been combined with their counterterms. The 
calculations for these contributions are performed with the help of FeynArts, FormCalc 
and LoopTools. The virtual photonic corrections (the diagrams are similar to the ones of 
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the leptonic decays, i.e. Figure are IR divergent. The contributions of the real photon 
bremsstrahlung are necessary for the cancellation of the IR divergences. We neglect the 
light quark masses, i.e. mg{q = u, d, s), except when they appear in the one- loop integrals. 
This gives a mass singularity Inm^. The virtual photonic corrections can be calculated 
analytically in the same way as in Section 14.2.21 . 

The diagrams for the real photon bremsstrahlung X2 ~^ XiQQI ^ire similar to the di- 
agrams in Figure 14. 5[ In analogy to Section I4.2.3[ the contributions of the real photon 
bremsstrahlung are also splitted into an IR-divergent soft part and an IR- finite hard part. 
For the light quark final states, we separate the hard photon bremsstrahlung into a collinear 
part and a non-collinear part in order to obtain stable numerical results. As presented in 
Section 14.2.31 the soft and collinear contributions are calculated analytically. The ana- 
lytical expressions for the singular part of the virtual photonic corrections, the soft and 
collinear contributions can be obtained by performing the replacements 

Q]^Q\, Inm; Inm, (4.57) 

in the corresponding expressions in Section 14. 2. 21 and 14. 2. 31 We treat the decays with heavy 
quark final states in the same way as t~t^ final states. The QED corrections are defined 
in the same way as in Section 14.2.41 since the photonic contributions are UV divergent and 
cannot be treated separately. 
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4.5 Numerical Results and Discussions 

In this section we present the numerical results both for a scenario where X2 can undergo 
two-body decays X2 ~^ ^^^^ ~^ and for a scenario where X2 o^^Y has three-body 

decays X2 ~^ Xi^~^~^ ■ The two-body decays are calculated in a complete and an approximate 
way, and the corresponding numerical results are compared. We discuss the total decay 
width of X2 and the branching ratios of its leptonic decays. The dilepton invariant-mass 
Mi+i- distribution is also presented and discussed, where the dilepton invariant mass Mi+i~ 
is defined as 

Mi+i- = ^{k+ + k-y. (4.58) 

As discussed in Section 14.2. 3[ the dilepton invariant mass Mg+e- is defined as coUinear- 
safe observable, i.e. we add the momentum of a coUinear photon to that of the emitting 
electron, since it is difficult to separate their energies at the LHC The energies of a muon 
and its collinear photon can be disentangled easily at the LHC, hence the dilepton invariant 
mass M^+^- is defined as non-collinear-safe observable, i.e. the momentum of a collinear 
photon is not added to that of its emitter muon. In this case the large logarithm In can 
not cancel in the distribution, so the mass effect can be seen in the dilepton invariant-mass 
distribution. In our calculations the selectrons and smuons have equal masses and the light 
lepton mass mi {I = e, /i) is neglected except when it appears in the one-loop integrals, so 
one will obtain identical distributions for M^+e- and M^+^- if both of them are defined 
as coUinear-safe observables. In order to see the differences of the two treatments (adding 
and not adding the momentum of a collinear photon to the emitting lepton), we also show 
the comparison of dilepton invariant mass M^+^- and M^+e- distributions. 

4.5.1 SPSla Parameter Set 

In the MSSM, soft-SUSY breaking is implemented by adding all possible soft terms to 
the Lagrangian instead of assuming a particular SUSY-breaking mechanism. All the pa- 
rameters in the soft-SUSY breaking Lagrangian fl2.52p are general matrices in flavor space 
and may be complex. This leads to more than a hundred unknown free parameters in 
the MSSM. Therefore, it is not practicable to scan over the entire parameter space. The 
"Snowmass Points and Slopes" (SPS), where several "benchmark scenarios" have been 
suggested [HIB], are meant to illustrate characteristic features of various scenarios of SUSY 
breaking. Among those, the so-called SPSla parameter set has been studied particularly 
widely. 

The SPSla parameter set is deflned in the framework of the mSUGRA scenario [50] . 
where the SUSY-breaking mechanism is supposed to be minimal supergravity. This sce- 
nario is characterized by four parameters and a sign, the scalar mass parameter Mq, the 
gaugino mass parameter M1/2, the scalar trilinear coupling Aq, the ratio of the Higgs vac- 
uum expectation values tan j3, and the sign of the supersymmetric Higgs mass parameter 
fi. These parameters are deflned in the DR scheme at the GUT scale Mqut- The SPSla 
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benchmark point is defined by setting [7j 

Mo = 100 GeV , Ml/2 = 250 GeV , Aq = -100 GeV , tan/3 = 10 /i > . (4.59) 

The low-energy parameters in the MSSM are obtained via renormahzation group running 
from the high-energy scale to the weak scale. This can be performed with various programs 

In Ref. the corresponding low-energy parameters for the SPSla benchmark point 
are obtained with IS A JET 7.58 [62], which read 

Mg = 595.2 GeV, tuao = 393.6 GeV, = 352.4 GeV, 

tan/3 =10, Ml = 99.1 GeV, M2 = 192.7 GeV, (4.60) 

where Mi and M2 are the gaugino mass parameters, Mg denotes the mass of the gluino, m^o 
is the mass of the neutral CP-odd Higgs boson A^. The soft SUSY-breaking parameters 
in the diagonal entries of the squark and slepton mass matrices have been chosen to be 
the same for the first and second generation. The off-diagonal entries have been neglected 
for the first two generations, i.e. there are no sfermion mixing. The soft SUSY-breaking 
parameters in the diagonal entries of the the squark and slepton mass matrices are 

Mq^^ = Mq2^ = 539.9 GeV, Mj^ = 519.5 GeV, M^^ = 521.7 GeV, 

M;^ = 196.6 GeV, Mj^ = 136.2 GeV , (4.61) 

where the index i in Mq-^ denotes the generation, u and d denote the up- and the down- 
squarks for the first two generations and I stands for the first and second generation slep- 
tons. The soft SUSY-breaking parameters in the diagonal entries of the squark and slepton 
mass matrices of the third generation have the following values, 

Mq3^ = 495.9 GeV, M^^ = 516.9 GeV, M^^ = 424.8 GeV, 

Mf^ = 195.8 GeV, M^^ = 133.6 GeV, (4.62) 

while the trilinear couplings of the third generation read 

At = -510.0 GeV, A = -772.7 GeV, A^ = -254.2 GeV. (4.63) 

The SPSla benchmark point gives rise to a particle spectrum where many states are 
accessible both at the LHC and at the ILC ^6j . The spectrum of supersymmetric particles 
at this benchmark point is shown in Figure 14.101 [U El ES]- Note in particular that the 
two-body decays X2 ~^ ifl^ Xi^^^^ kinematically allowed. No other two-body 
decay mode for X2 is open. Moreover, the sleptons have the same masses for the first two 
generations, i.e. m^;^ = m£>^, mg^ = m^^, mg^ = m^^. 

Another set of parameter is called SPSla' which is proposed in the SUSY Parameter 
Analysis (SPA) project [66]. Its parameters are close to the snowmass point SPSla except 
for a small change of the scalar mass parameter and the trilinear coupling to comply with 
the measured dark matter density. The low-energy parameters for the SPSla benchmark 
point which have been shown in fl4.60[ 14.611 14.611 14.631) are used as the input parameters 
in calculating the decays of X2- 
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Figure 4.10: The spectrum corresponding to SPSla benchmark point 



4.5.2 Numerical Results for the SPSla Parameter Set 

The SPSla parameter set has been discussed in Section 14.5. H where the mass spectrum 
of the MSSM particles was given in Figure I4.10[ Since the two-body decays X2 ~^ 
are kinematically allowed, the decays X2 ~* Xi^~^^ ^ire calculated in both complete and 
approximate way as discussed in Section l^y2] and lT3l The numerical results of both methods 
are compared. 

We calculate the QED corrections of the decay X2 ~^ Xi^^e^ for the SPSla parameter 
set and show the contributions F + Fcou, Fnon-cou and Fqed as function of AE and A9 
in Figure 14. Ill From these figures one obtains that the QED corrections do not depend 
on these cut-off parameters when they are very small. This is in accordance with the 
discussions in Section 14.2.31 that the cut-off parameters should be small enough so that 
the soft and collinear contributions can be treated approximately as in (I4.34p and (14.401) . 
The instabilities of QED corrections can be seen when the values of the cut-off parameters 
become too small. As a result, we use A6 = 1° and AE = O.OOl-^/i in our calculations. 
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Figure 4.11: The dependence of F + Fcou, Fnon-cou and Fqed on the cutoff parameters 
AE and A6, where we have used AO = 1° and AE = 0.001-\/i in the left and right plot, 
respectively. Here ^/s = m^o denotes the center-of-mass energy. The results are from the 
decay X2 ~^ Xi^'e^ for the SPSla parameter set. 



Dilepton Invariant Mass Mi+i- Distribution 

The dilepton invariant mass M^+^- distribution from the complete calculation is shown in 
Figure 14.121 In the left frame we show not only the tree- and one-loop-level predictions, 
but also the separate QED and "remainder" corrections, see fl4.48p . We see that the non- 
QED contributions are positive and quite large everywhere, whereas the QED contribution 
is large and negative near the endpoint of the distribution, but small elsewhere. In full 
three-body kinematics this endpoint is simply given by 



e+e l3— body 



^2 



80.4GeV , 



(4.64) 



where the numerical value holds for the SPSla scenario. Moreover, in this scenario X2 
decays are dominated by contributions with on-shell li in the intermediate state. The 
endpoint for this two-body configuration is given by 



-body = '^X§' 



2 



1 - 



1 - 



2 



76.8 GeV 



(4.65) 



Note that this is only 3.6 GeV below the endpoint of the three-body decays. At tree 
level, the Mg+g- distribution peaks at the region which is a little below the endpoint of 
the two-body contribution. The right panel in Figure I4.12[ which shows a blow-up of 
the peak region, shows that the shape of Mg+e- distribution altered once higher order 
corrections are included, i.e. the peak of this distribution is moved downwards by about 4 
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Figure 4.12: The dilepton invariant mass Mg+g- distribution for the SPSla parameter set 
(complete calculation). 



GeV. This is almost entirely due to contributions where a hard photon is emitted, which 
takes away energy from the e^e~ system. This change of the shape of the invariant-mass 
distribution near the endpoint is important, since in (simulated) experiments one needs 
a fitting function describing this distribution in order to determine the location of the 
endpoint [67] . 

In Figure 14.131 we compare the numerical results of the complete calculation and the 
single-pole approximation at tree (left) and one-loop level (right). At tree level the Me+e- 
distribution computed in the single-pole approximation has an exactly triangular shape, 
with a sharp edge at the endpoint fl4.65p . This edge is smeared out a bit in the com- 
plete tree-level calculation, which includes complex slepton masses in the propagators as 
fl4.7p . As noted above, this edge is also softened considerably once hard photon emission 
is included. The single-pole approximation therefore works even better in the one-loop 
calculation. However, this excellent agreement even for the differential decay width is par- 
tially accidental. The agreement would become somewhat worse if the endpoints in two- 
and three-body kinematics were further apart. This would happen if the mass of l\ was 
close to the mass of either X2 oi' Xi? since then one of the two square roots in fl4.65p would 
become small. 

The comparison of the dilepton invariant mass M^+^- and Me+e- distributions are 
shown in Figure 14.141 In the upper frames we show the dilepton invariant mass Mi+i- {I = 
e, /i) distribution both at tree and one-loop level. Since the selectrons and smuons have 
equal masses and the light lepton mass rrii (/ = e, /i) is neglected except when it appears in 
the one-loop integrals, their distributions are identical at tree level and different at one-loop 
level due to the different treatment of the collinear-photon radiation. From these figures 
one obtains that at one-loop level the mass effect is larger near the endpoint than in other 
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Figure 4.13: The dilepton invariant mass M^+e- distribution for the SPSla parameter set 
(comparison of the complete calculation with the approximate calculation). 



regions and the peak of the M^+^- distribution is shifted to lower invariant-mass values in 
comparison with the Mg+g- distribution. We also show the relative one-loop corrections 
in the lower frames in Figure 14.141 The relative one-loop corrections from the fi~^fi~ final 
state is smaller than that of the e~^e~ final state in the upper invariant- mass region, while it 
is larger in the lower invariant- mass region. The main reason is that we add the momenta 
of coUinear photons to that of emitting electrons, but we do not do this for the coUinear- 
photon radiation from muons. Hence the invariant mass M^+^- is reduced in comparison 
with Mg+e-. This leads to the shifting of events from the upper invariant-mass region to 
the lower invariant-mass region. 

The corresponding results for the dilepon invariant mass Mt-+t— distribution are shown 
in Figures. 14.151 and 14.161 The left plot in Figure 14.151 shows not only the tree- and 
one-loop- level results, but also the QED and "remainder" corrections. The definitions 
were presented in Section 14.2.41 The endpoint region of the Mt-+t— distribution from the 
complete calculation is shown in the right panel of Figure 14.151 We see that the peak of 
the distribution is shifted downwards by about 2 GeV once higher-order corrections are 
included. A shift of this magnitude may be significant, even though the r+r" invariant 
mass is in general difficult to measure accurately, due to the presence of (anti-)neutrinos 
in the r decay products, which escape detection. In Figure 14.161 predictions from the 
complete calculation are compared to those from the single-pole approximation. In this 
case we find almost perfect agreement even in the endpoint region, both at tree level and 
after including one-loop corrections. The reason is that for the SPSla parameter set, 
happens to be very close to ^m^om^o. Performing the replacement mg^ — >■ ruf^ in (14.651) . 
one finds that the endpoint of the t^t~ distribution in two- and three-body kinematics 
practically coincide. 
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Figure 4.14: The comparison of the dilepton invariant mass M^+^- and Mg+g- distribution 
for the SPSla parameter set. 
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Figure 4.15: The dilepton invariant mass Mt-+^- distribution for the SPSla parameter set 
(complete calculation). 




Figure 4.16: The dilepton invariant mass M,-+^- distribution for the SPSla parameter set 
(comparison of the complete calculation with the approximate calculation). 
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Total Decay Width of X2 the Branching Ratios of the Decays X2 ~^ Xi^~^^ 

The partial widths of the different X2 decay modes and the branching ratios of its visible 
leptonic decays are listed in Table 14.11 where the numbers in the parentheses are obtained 
from the approximate calculations. We find ^iT{x2 ~^ Xi^i^i) J2i^iX2 ~^ Xi^~^^~)- 
This is not surprising, since the charged lepton final state is accessible via on-shell /i 
intermediate state, whereas for the neutrino final state all exchanged particles are off shell. 
Since squark masses are near 500 GeV in SPSla scenario, hadronic final states contribute 
even less than neutrinos do. 



decay mode 


tree-level width(MeV), Br 


one loop-level width(MeV), Br 


e'e+Xi 


1.123 (1.122), 5.9% 


1.297 (1.294), 6.7% 




1.123 (1.122), 5.9% 


1.297 (1.294), 6.7% 




16.870 (16.933), 88.0% 


16.595 (16.646), 86.2% 


l^ei^eXl 


0.012 


0.012 


J^fiJ^t^Xl 


0.012 


0.012 


PrT^rXl 


0.013 


0.013 


QQXl (Qi^t) 


0.015 


0.015 


total width 


19.168 


19.241 



Table 4.1: Partial widths of different X2 decay modes and the branching ratios of its visible 
decays for the SPSla parameter set. The numbers in parentheses give the corresponding 
partial widths calculated in the single pole approximation. 

From the results in Table 14.11 one concludes: 

• The main decay mode of X2 is X2 ~^ '^"'^'^Xi- Its branching ratio is about 88.0% at 
tree-level, 86.2% at one-loop level. This mode dominates partly because of the lower 
mass of fi as compared to ei (133.0 GeV vs 142.7 GeV). Even more important is 
that li is a pure SU{2) singlet for I = e, fi, since we neglect terms oc mi in the mass 
matrices of these sleptons. In contrast, — tji mixing is quite significant, leading to 
a sizable SU{2) doublet component of fi. Therefore X2 decays into (real or virtual) 
h can only proceed through its small U{1)y gaugino (bino) component for / = e, /i, 
while the large SU{2) gaugino (neutral wino) component also contributes for / = r. 

• The total X2 decay width is enhanced by 0.4% when one- loop corrections are included. 
Such modest corrections are typical in the absence of large enhancement factors (e.g., 
large logarithms). 

• One-loop corrections enhance the partial width and the branching ratio of X2 ~^ 
l"l^Xi — ^'A*) decays by 15.5% and 13.6%, respectively. 

• The single pole approximation reproduces the integrated partial widths to about 
0.3% accuracy. This agreement is even better than in the Mi+i- distribution shown 
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in Figures 14.131 and 14.161 From (14.81) and the discussions of the large logarithm 
In (|r[^^V"^[i) in Section 1122] one might expect better agreement for the integrated 
partial width than for the kinematical distributions. 

4.5.3 Numerical Results for the Pure Three-body Decays 

We also investigated the effect of higher-order corrections on leptonic X2 decays for a 
scenario where X2 does not have any two-body decay modes. To that end we again use the 
SPSla parameter set, except that the soft SUSY-breaking parameters in the slepton mass 
matrix are set to 

rrii^ = 230GeV , m^^ = 183 GeV, I = e, fx, r . (4.66) 

The masses of the relevant neutralinos and sleptons in this modified SPSla parameter set 
are listed in Table 14.21 where one finds that X2 has to undergo a pure three-body decay. 
Therefore we do not have to introduce complex slepton masses in the one-loop functions. 



particle 






ei (/ii) 


62 (1^2) 


n 


T2 


ui {I = e,fi,T) 


mass (GeV) 


176.6 


96.2 


187.9 


234.9 


182.3 


239.2 


221.0 



Table 4.2: Masses of the relevant neutralinos and sleptons for the modified SPSla 
parameter set 

The dilepton invariant mass M^+e- and Mt-+^- distributions are shown in Figures [4. 171 
and 14.181 respectively. At tree level the Mg+e- distribution shows a small peak near its 
upper endpoint from the exchange of nearly on-shell Z bosons. Since the QED and non- 
QED corrections are very small and negative in this region, this peak is less pronounced 
once one-loop corrections are included. This is of some significance, since the shape of 
this distribution can now be used to infer the strengths of various contributing diagrams, 
which in turn provides information on slepton masses and neutralino mixing pUl [55] . Since 
f exchange is much enhanced relative to e exchange, one cannot see any contributions of 
Z— exchange even at tree-level from the M^+t— distribution. Moreover we can observe 
that the invariant mass Me+e- and M-^+r- distributions have a rather sharp edge at their 
endpoints. These edges are again softened by real photon emission, but remain quite 
distinct. This should facilitate the experimental determination of the endpoint, and hence 
the measurement of m^o — m^o . 

A2 Al 

We compare the dilepton invariant mass M^+^- and Mg+e- distributions in Figure 
I4.19[ where the tree- and one-loop-level results, the blow-up of the endpoint region and the 
relative one-loop corrections are shown. From these figures one obtains that the shapes 
of the M^+^- and Mg+e- distributions are identical at tree level and different at one-loop 
level due to the different treatment of collinear-photon radiations. In contrast to the 
numerical results from the SPSla parameter set (see Figure 14.141) . the mass effect is small 
in Figure 14.191 but it is still distinct, especially in the relative one- loop corrections. In 
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Figure 4.17: The dilepton invariant mass Mg 
three-body decay. 
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Figure 4.18: The dilepton invariant mass M^+t— distribution in the case of a genuine 
three-body decay. 
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the calculations for the invariant-mass distribution, the momentum of a collinear photon 
is added to that of the emitting electron, but it is not added to that of the emitting muon. 
Hence the invariant mass M^+^- is reduced in comparison with Me+e-. It leads to the 
shifting of events from the upper invariant-mass region to the lower invariant-mass region. 
This effect can be seen in the lower frames in Figure I4.19[ i.e. in the lower invariant-mass 
region the relative one-loop corrections of the final state is larger than that of e^e~ 

final state, while the inverse relation holds in the upper invariant-mass region. 
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Figure 4.19: The comparison of the dilepton invariant mass M^+^- 
in the case of a genuine three-body decay. 



and Mp+p- distribution 
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The partial widths of different X2 decay modes and the branching ratios of its visible 
decays are shown in Table 14.31 The t~^t~Xi final state is still the largest decay mode 



decay mode 


tree-level width(keV), Br 


Hoop-level width(keV), Br 




1.270, 4.4% 


1.451, 5.0% 




1.270, 4.4% 


1.451, 5.0% 


r-r+x\ 


7.209, 25.1% 


7.383, 25.4% 


^el'eX\ 


1.273 


1.355 




1.273 


1.355 


Ur^rXl 


1.273 


1.354 


UUXi 


2.480 


2.386 


ddxi 


3.330 


3.298 


ccxi 


2.475 


2.378 


ssxi 


3.330 


3.298 


bbxi 


3.595 


3.405 


total width 


28.778 


29.114 



Table 4.3: The decay width of different X2 decay modes and the branching ratios of its 
visible leptonic decays in the modified SPSla scenario. 

of X2 (25.1% at tree level, 25.4% at one-loop level), but it no longer dominates. The 
hadronic final states have very large partial decay widths and branching ratios: r^^J^^j.jj^i^ = 
15.210keV (52.9%), Tl'^^ic = 14.765keV (50.7%), though the squark masses are much 
heavier than the slepton masses. The part reason is that the Z-exchange diagrams give 
larger contributions to hadronic final states than to leptonic ones. Moreover the interference 
between Z and squark exchanges is large and positive for the hadronic final states, while 
the interference between Z and slepton exchanges is also large but negative for the leptonic 
final states. This is the main reason why the hadronic decays of X2 obtain so large branching 
ratios. 

Note that exchange of the SU{2) doublet sleptons now dominates for I = e, fi. This 
dominance of exchange also explains why the e~^e~Xi and z/eZ/gX? final states now have 
quite similar partial widths. We have assumed three exactly degenerate sneutrinos here, 
unlike in the original SPSla scenario, where 9^- is slightly lighter than Uf,. In the modified 
scenario a tiny difference between the one-loop partial widths for VrV-rX^ and Uei^eXi final 
states nevertheless results from one-loop corrections involving the r mass or Yukawa cou- 
pling (e.g. from the u and u two-point functions). The total X2 decay width is increased 
by 1.2% when one- loop corrections are considered. The partial width of the X2 decay into 
electron and muon pairs is now enhanced by about 14.3%, leading to an increase of the 
corresponding branching ratios by 13.6% at one-loop level. 



Chapter 5 
Conclusions 



In the MSSM with conserved R-parity, the decays of the next-to-hghtest neutrahno X2 i^^o 
the LSP Xi and two fermions are always involved in the decay chains of supersymmetric 
particles. Moreover, X2 '^^^ of lightest visible supersymmetric particles that can 
be produced directly at future e+e" colliders and plays a prominent role in the analysis 
of cascade decays of gluinos and squarks at the LHC. An accurate understanding of its 
decays is therefore of considerable importance. The leptonic decays of X2 particularly 
interesting since the endpoint of the dilepton invariant mass distribution can be used to 
reconstruct the mass differences of the supersymmetric particles. In this thesis, we have 
investigated X2 leptonic decays X2 ~^ Xi^~^~^ one-loop level. 

For the cases where the intermediate charged sleptons h can be on shell, these decays 
were calculated both completely and in a single-pole approximation at one-loop level. In the 
complete calculation one has to employ complex slepton masses in the relevant propagators 
and one-loop integrals. The single-pole approximation in this case is performed in the way 
that the X2 decays are treated as a sequence of two two-body decays. For the numerical 
evaluation we use the SPSla parameter set. We compare the results from the complete 
and approximate calculations and find that this approximation reproduces the integrated 
partial widths to better than 0.5% accuracy even after one-loop corrections are included. 
From these calculations one obtains a rather small one-loop correction to the total X2 
decay width, but the branching ratios for the electron and muon final states are increased 
by about 13.6% at one-loop level. 

The dilepton invariant mass Mi+i- distributions were also studied. The shape of these 
distributions is found to be altered by real photon emission contributions, i.e. its peak 
is shifted by several GeV below the endpoint. This is very important since the shape of 
the distribution near the endpoint should be known if the endpoint is to be determined 
accurately from real data. In our calculation we define coUinear photons as being emitted 
at an angle < 1° relative to the emitting lepton. Since the selectrons and smuons 
have equal masses and the light lepton mass m; {I — e, /i) is neglected except when it 
appears in the one-loop integrals, one will obtain identical distributions for M^+e- and 
M^+^- if the momentum of a collinear photon is added to that of the emitting lepton. 
The actual effect of the coUinear-photon radiation depends on details of the measurement 
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apparatus, and therefore has to be calculated anew for each experiment. We have focused 
on the LHC experiment in our calculation. At the LHC the electron energy is determined 
calorimetrically. In this case a coUinear photon would hit the same cell of the calorimeter as 
the electron, so the two energies cannot be disentangled. Hence we add the momentum of 
a collinear photon to the one of the emitting electron in our calculation. Since muons pass 
through the calorimeter, where the photons are detected, and measured forther outside in 
the muon detector at the LHC, the momentum of a collinear photon is not added to the 
one of its emitter muon in our calculation. In this case the mass effect can be seen in the 
dilepton invariant-mass distribution. We find that the peak of Mg+g- distribution is moved 
downwards by about 4 GeV once the one-loop corrections are added. In contrast to the 
Mg+e- distribution, the peak of the M^+^- distribution is a little shifted to lower invariant- 
mass values at one-loop level. This is due to the different treatment of the coUinear-photon 
radiation. 

We have also analyzed a scenario with increased slepton masses, so that X2 can only 
undergo genuine three-body decays. We find that the total X2 decay width is enhanced 
by 1.2% when one-loop corrections are considered, while the branching ratios of X2 decay 
into electron and muon pairs are enhanced by about 13.6% at the one-loop level. One 
also finds that the shape of dilepton invariant mass Mi+i- distributions are also affected 
by the real photon emission. Moreover, these distributions have a rather complicated 
shape, showing the contributions from Z exchange near the upper endpoints. In this 
case the shape of the distribution away from the endpoint also carries information about 
slepton masses and neutralino mixing angles. Fitting tree-level distributions to real data 
might therefore give wrong results for these physical parameters. In this context a careful 
analysis of the collinear- photon radiation is also important, since differences in the energy 
measurements of electrons and muons could lead to spurious differences of fitted selectron 
and smuon masses. Here the collinear-photon radiations for electrons and muons are 
treated as discussed beforehand. One finds that the one-loop shapes of the Mg+g- and 
M^+n- distributions are different, though the selectrons and smuons have equal masses in 
our calculations. 



Appendix A 

Notations and SM Parameters 



In this thesis we adopt standard relativistic units, i.e. h — c — 1. A general covariant 
four- vector is denoted by 

an = (ao, ai, 02, 03) = (ao, a) (A.l) 
and a contravariant four- vector is 

= (a°, a^ a=*) = (oo, -a) . (A.2) 
They are connected by the metric tensor 

9^" = y^. = diag (1,-1, -1,-1) (A.3) 

via the relations 

a" = g'^'^a,. (A.4) 
The product of the four- vectors are defined as 

— * 

ab = a^bn — Oofco — a, ■ b . (A. 5) 
The four-gradients and are defined 



d 


i- 




\dt 


d 


(- 




\dt 



We also use the compact "Feynman slash" notation 

(i^ra,, (A.7) 

where 7^ are Dirac matrices. 
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A.l Pauli and Dirac Matrices 

The pauli matrices are defined as 

i_/Ol\ 2_/0-i\ 3 f 1 



They satisfy the commutator relation 

[a\a^ = 2i£^^V, i,j,k^ 1,2,3. (A.9) 
The totally antisymmetric tensors in three dimensions are defined as 

{+1 , for even permutations of 123 
—1 , for odd permutations of 123 (A. 10) 

, otherwise . 

One can arrange the Pauli matrices as 

= (a^a) = (a^a^a^a3) , (A.ll) 
= {a', -a) , (A.12) 

where '^^ ~ ^ q i ^ ' Anti-symmetric matrices a^^ and a^^^ are defined by 

a''" = ^ {a^a" - a^a^") , (A.13) 

= 1 (^^^cr'^ _ ^V'^) . (A.14) 

The Dirac 7-matrices are defined via the anticommutation relations 

{7^7'^} = 2^^^ (A.15) 

A fifth 7-matrix is defined by 

7^=75 = i7V7V- (A. 16) 

From these definitions one can easily obtain the following properties for the 7-matrices, 

{75,7/'} = 0, (7^)^ = 0. (A.17) 
In the chiral or Weyl representation the explicit expressions for the Dirac 7-matrices are 













j 




on 



7''= \ ,7^= n ; ■ (A-18) 



The left- and right-handed operators are defined by 



c^L = ^ (1 - 75) , = ^ (1 + 75) ■ (A.19) 
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A. 2 Spinors 

The components of the two-component (Weyl) spinor are Grassmann numbers, i.e. 

{Xa,Xp} = = =0, 

{Ca,Cp} = } = {ed,f } = 0, (A.20) 

and they also have anticommutation relations with other Grassmann numbers. Here the 
indices a{a) — 1,2 and = 1,2. The scalar product of two-component spinors x 
^ is defined as 

— X 1 
X^ = Xa^ 5 

Xa^^ ^ xV^r. (A.21) 
A four-component (Dirac) spinor ^ in the Weyl representation can be constructed via 

* = (|d ) , (A.22) 

where ^q, and are Weyl spinors. The Dirac-adjoint spinor ^ is expressed as 

* = *^7o = ( X" ) ■ (A.23) 
The charge conjugation of the Dirac spinor ^ is defined via 

= C^^ ^(^^2^ , (A.24) 
where the charge conjugation matrix C is expressed as 

A Dirac spinor ^ is also a Major ana spinor if the relation ^ = ^'^ is satisfied. Hence a 
Major ana spinor A can be written as 

A - (I ) . (A.26) 
The left- and right-handed components of a Dirac spinor can be written as 

*L = u;l^=(^ ) , (A.27) 

^u;r^=( ) . (A.28) 
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Some useful relations between the four- and two-component spinors are 



= X16 + 6X2, (A.29) 

^17^*2 = 6a''6-X2^'^Xi, (A.30) 

^i7'*2 = -X16 + X26, (A.31) 

^17^*2 = -Ci^'^^ - X2a'^Xi , (A.32) 

*i7''9^*2 = Xi<7'^9^X2 + 6^^9^6, (A.33) 

^ic^L*2 = X16, (A.34) 

= 6X2, (A.35) 

*i7'^a;i*2 = Ci^'^^, (A.36) 

*i7^a;fl^'2 = -X2^''Xi, (A.37) 

*i7'^a;z,9^*2 = ^^'^9^6 , (A.38) 

*i7'^a;«9^*2 = Xi<7'^9^X2 ■ (A.39) 



A. 3 SM Parameters 

For the numerical evaluation, the following values of the SM parameters are used: 

me = 0.510999MeV, = 105.6584MeV , = 1.777GeV , 

ruu = 53.8MeV, mc = 1.5GeV, mt = 175GeV, 

rud = 53.8MeV , m, = ISOMeV , = 4.7GeV , 

mw = 80.45GeV, = 91.1875GeV , 

a = 1/137.0359895, = 1.1663910 x 10-5GeV-^ 



Appendix B 
One-loop Integrals 



B.l Definition of the One- loop Integrals 

We define the one- loop integrals in the same notation as in LoopTools. As discussed in 
Chapter [3], dimensional reduction is used, where only the momenta are calculated in D 
dimensions, while the fields and the Dirac algebra are kept 4-dimensional. The definition 
for the scalar one-loop integrals are shown in the following, their tensor integrals are defined 
by adding the momenta g^, g^gj, ■ ■ ■ to the numerator. 




(27r/x)^-^ r 

in'^ J [g2 - mf] [(g + - m|] [(g + pi +^2)^ - "^1] ' 
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Do {pI^pIpIpI (Pi +P2)^ {p2+P3f,mj,mlml,ml) 



m'^ J [g2 _ ml] [{q + pi)^ - m|] [{q + pi + P2Y - mf\ [{q +P1+P2+ Pzf - "t-I] ' 



B.2 Scalar One- loop Integrals 

The general formula for the scalar one-, two-, three- and four-point functions were derived 
in [5lj . Here we not only outline their general formula but also give the explicit expressions 
in some special cases. 



Scalar One-point Function 

The scalar one-point function can be written as 



Aoim") = ( A - In ( --^ 1 + 1 



where the UV-divergent part A is defined via 

2 



4-D 



7^; -|- In 47r 



with 7e is Euler's constant, and /x is renormalization scale. 



fB.r 



(B.2) 



Scalar Two-point Function 

The scalar two-point function can be written as 



Bq{p^ , ml, ml) = A — / dx\n 

Jo 



X p — x{p — ml + mi) + ml — ie 
7 . 



(B.3) 
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where e is a infinitesimal real number. Below are some special cases, 



Bo{mf,mf,0) = A - In 



mr 



+ 2 



,mt 



5o(0,m ,0) = A-ln(^) + l, 



+ 2, {p^ > m; 



So(p^m,^m|^) = A-ln 
5o(/,0,m?) = A-ln | ^ ) +2, (p^ ~ ) . 



mi 



(B.4) 
(B.5) 
(B.6) 

(B.7) 



Scalar Three-point Function 

The scalar three-point function can be expressed as 



Co {pI, pI, {Pi + P2)^ ml, ml, ml) 
where 



6{1 - Xi- X2- Xs) 



dx I dx^ — r — 7" 

{Xi+ X2+ X3)g{Xi,X2,Xs) 



(B.8) 



g{xi,X2,Xs) = {mlxl + mlxl + mlxl){xi + X2 + xs) - 
pIxiX2 - plx2X'i - (pi + ^2)^2:1X3 - ie . 



(B.9) 



A special case can be expressed as 
Co {nn?i, mf, (pi + p2)^ 0, mf, mf) 



1 



1 



In 



mf 



-{Pi +P2y - 



In ( 4 1 + 



mf 



m 



-In', . 

2 V.-bi+Psj^-^e 



TT 



(B.IO) 



where A is the photon mass regulator, and we have assumed (pi + ^2)^ ^ m\ . This can 
also be found in Ref. [55j. 

In the calculations for ^ decays one has to calculate the scalar three-point function 
Cq = Co{mf,m'io, (pi +p2Y,0,mf,mf) {{pi + P2Y ^ mf) analytically. It can be obtained 



via the calculation of fIB.Sp . Below are the analytical expressions in different cases. 
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For the general X2 three-body decays, i.e. (pi + P2Y 7^ mf , it can be expressed as 

is 

follows, 

1 [ 1, ( ^ 1, ( ^-'^-'^ ] , 

° (Pi + - m'io \ \{Pi+ - rn'io - «e y ^"^1^ - {pi + P2Y - y 



2 / mj - {pi + - ie\ ( {-mlQ - ie){m} - {pi + P2Y - ie) 



+ ^ - Lis I '° o ^— I + 



6 I mf — m?o — j \ (— ""^-o + (Pi + ^2)^ — "^e 



-Li2 



-m-o — ie 



{Pi+P2y - - ie 



(B.ll) 



The dilogarithm Li2{x) is defined as 



Li2{x) = - f dt^^^-^ . (B.12) 
Jo i 

When the sleptons Ig can be on shell, i.e. {pi + P2Y is close to m?- , using complex 
masses mf- — iTr m,- , one obtains 

Is Is 's 



{Pi + P2Y - mlo I 2 \mf - m|o - «e / 6 



X2 K \ Is X2 

(pi + P2)^ - m'lo ~ '^^ ) y^} ~ (^'1 + P'^y ~ '''^is'^is ~ "^^ 

m} \ f {pi+ P2Y - - ie' 

In I ■ In ' 



(B.13) 



- {Pi + ^2)^ - ^r^- mf - ie j \ {pi+ P2Y 



Li 2 



-mln — ie 

xi 



{Pi +P2y - - ie 

X2 



Here we have used Li2{0) = 0. This formula is for the complete calculation. 

In this case one can also treat X2 decays as production and decay of Is, i.e. {pi +^2)^ 
m^- , the function Cq is IR divergent. The explicit expressions are 



mf — mlo I \ m"- — mtg — le j Xmim,- I \ m} — mt^ — le 



Is x'i K \ Is x'i / y-i, 

.2 



■\^^\^r^2 — -pM^A — -1 ^ (B.14) 

2 \ mf — mio ~ y ~ ~ "^^ 



which is also given in Ref. [55] . 
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Scalar Four-point Function 

The scalar four-point function DQ^mf, m^^rnio^mj, {pi +^2)^, {p2 +^3)^, 0, m?, mf , m?) is 
necessary for the calculation of the three-body decays of X2- Its analytical expressions can 
be obtained from Ref. [55] , 



1 1,0/ i^f^i^i 



iP2 +Piy{m} - {pi +P2] 



In' 



mf- — mln - It 



+ ln" 



nil mi 



mf- — m\ — "Je 

Is Xl 



+ 



2 In 

3 



mf 



In 



mr X 



-(P2 + -i^J \m} - {pi + p2y - inif Tr - it 



+ 



(m?„ — m^- + ze)(m-o — + ^e) 



(B.15) 



Here we focused on the case where {pi + ^2)^ is close to m?- and used the complex masses 



mj- m,- — iVr mr . 

Is Is '-^ 



Photonic Part of the Fermion (Sfermion) Field Renormalization 
Constants 

The photonic part of the fermion self-energies are given by 



photonic 
photonic 



-^6,,Q}[2B,{p\m%0) + l], 



a 
'2^ 



photonic AtV 

6,,Q}mf,[2Bo{p\ml,0)-l]. 



(B.16) 



The fermion field-renormalization constants have been presented in f l3.25bp and fl3.25cl) in 
Chapter [31 In order to calculate them in the case of light fermions we need the relations 



dBo{p\mlO) 



o 9 



dBi{p'^,mj,0) 



1 A 
2 In 



m 



1) 



dp^ 



o 9 



2(M^)-A-3) 



m 



(B.17) 



96 APPENDIX B. ONE-LOOP INTEGRALS 

Here we have used the general relations shown as follows (see also 

2 2 1 

5i(p^mi,m2) = -^—-^{Bo{p'^,mi,m2)-Bo{0,mi,m2))--Bo{p'^,mi,m2), 

2p^ 2 

dp^ ~ — 2p^ — ^^^P ' ~ '^^^ ^ 

^2_^2_p2 dBo{p\m,,m2) 

dBo(p'^,mi,m2) m^ — m^ , m2 mim2 (\ 

— m 1 — r m r 



dp^ p4 p4 

p^ \ — 1 



3^ l + -^lnr), (B.18) 



where r and ^ are determined from the equation 

^2^^ 2 — + i = + + (B.19) 

mxm2 r 

The photonic part of the field renormalization constants for the light fermions can be 
expressed as 

^^'''Iphotonic = ^^•''Iphotonic = (in - 4 lu - A - 4^ . (B.20) 
The photonic part of the sfermion field renormalization constants can been written as 



^Zi\ , , . = -Re 



2^ p'^=m?. 

, (B.21) 

photonic 



fa I photonic 3p^ 

where the sfermion self-energies are expressed as 

S/Jp^) = -^gJ/(So(p^O,m|)-^?l(p^O,m|)) . (B.22) 



Using the relations in (IB.ISP one obtains the expression for The IR-singular 

part reads 

'^'^/ii I photonic yj-^/^''^ I 777, J ^ ' (B.23) 
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Multi-Channel Monte Carlo Method 



C.l Principles of the Monte Carlo Method 

The Monte Carlo method (see Ref. [02]) is a way to calculate the integrals with a large 

number of integration variables. The integration variables x are mapped to a set of random 
numbers r via 

x = h{r), 0<r<l. (C.l) 

An integral can be written as 

fi^)dx= I'^^dr, (C.2) 
Jo 9{h{r)) 



where g{h{r)) is density, it is defined as 

1 dh{r) 



g{h{r)) dr 



(C.3) 



This integral can be approximated by sampling the integrand times and taking the 
average, 

N £^ g{h{n)) ■ ^ ■ ' 

The integration error is defined by 



P -P 

(C.5) 



N 

When the integrand / varies strongly in the phase space, the efficiency of the Monte 
Carlo method is improved by Importance Sampling, i.e. more events are sampled in the 
important region where / becomes large. This is implemented by choosing the variables x, 
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the mappings between these variables and tlie random numbers r in such a way that the 
resulting f/g is much smoother than /. 

In practice, we choose the variables in such a way that the Lorentz invariants corre- 
sponding to the propagators are included. The mapping is chosen such that the density 
g behaves in a similar way as the propagator. Below are the mappings belonging to the 
different propagator types [58] . 

• Propagator with vanishing width ^— 

The explicit expressions for the mapping x = h{r) and the resulting density are 
written as follows, 

/i(r,m^ z/,Xmm,a;max) = ['^(a^max - m'^Y'" + (1 - r)(xmm - m^y^"]^ 
for < X and 7^ 1, 



/l(r,m^, Z/,Xmm,a;max) = -[^("^^ - Xmax)^ " + {1 - r){m^ - XminY ""Y-" ^ VT? , 
gx{x,'m , Z^, Xinin; a;niax) — 77 o 7 o o \v (^•'^) 

[{m^ - Xmax)^ - [m^ - x^n^Y '^J (m^ - xf 

for m? > X and z/ 7^ 1 and 



h{r, z/, Xmin, a;max) = exp [r ln(a;max - rn^) - (1 - r) ln(a;min - rn^)] + 

gx{x,m , Z/, Xmin; a^max) = Tj z ;7^ ; Z oTTT oT (C.^ 

[ln(xinax - m^) - ln(Xmin " m^)\ [x - m^) 



for u = 1. 

• Breit- Winger propagator ^_^2^-^-p : 
The variable x is mapped to 

/i(r, - imT, 2, Xmin, a^max) = tan[?/i + (?/2 - yi)r] + , (C.9) 

the resulting density is 

2 niT 

gx{x^1Tl ZTTiF, 2, Xmin; a^max) 7 777 :7v5 ; 9-^91 5 (C.IO) 

(1/2 — 1/1) [(a^ — ^ ) + m^T^J 

where 

yi/2 = arctan( '"'"/'"'^^ ) . (C.ll) 

mi 

The parameter z/ can be tuned to optimize the Monte Carlo integration and should be 
chosen u > 1. Other variables, i.e. the polar and azimuthal angels 6 and 0, are generated 
as following, 

cos^ = 2r-l, = 27rr. (C.12) 
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C.2 Kinematics 



For a mutiparticle process, i.e. a particle with momentum ki ( mass mi) decays into n 
particles with momenta ki (mass rrii), the phase-space element is given by 



2F 



1=2 



iO 



S\k, 



n+1 

E 

i=2 



(C.13) 



This process can be described by Sri — 4 independent variables. In order to obtain the 
kinematics we treat the multiparticle process as taking place via cascade decays, where 
the intermediate states are unstable particles which then decay to others and eventually 
form the final states particles. Hence a mutiparticle process can be composed by isotropic 
particle decays, which are described as follows. 

One particle with momentum ki decays into two particles with momenta k2 and k^, 
masses m2 and m^. The polar angle (p aiid azimuthal angle 6 in the rest frame of the 
decaying particle are chosen to be the suitable integration variables. The phase-space 
element is defined as 



d^{ki, ml, ml) 



d^ko d^k-i 



2k2o 2fc 



5^^\ki-k2 



30 



\^l\klmlml) 
8kl 



2n 



1 



(C.14) 



d cos 6 , 



A is defined as follows, 



X{x,y,z) 



+ + — 2xy — 2yz — 2xz 



(C.15) 



Using the Monte Carlo method, the angles have to be mapped to the random numbers as 
flC.12p . hence the density can be written as 



n 2 2 2 

9d{ki,m2,m^ 



2kl 



(C.16) 



Since the laboratory frame usually does not coincide with the rest frame of the decaying 
particle, the Lorentz transformation is introduced [TD]. The Lorentz transformation of 
momentum k into the rest frame of the particle with momentum p is defined by 



where 7 



po 



,7/5 



b1 



and m 



k' = B{^,P)k, (C.17) 
'p2- the explicit formula for B{'~f,i3) can be written as 



( 1 7/3 \ 

10 

10 

V7/? 7 / 



(C.18) 
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The inverse Lorentz transformations is defined by replacing p by —p, 



(C.19) 



The orientation of the coordinate system can be arbitrarily chosen because the decay 
is isotropic. The momentum of the outgoing particle can be written as 



k = R{(j),9)B{-f,-p)k' 



(C.20) 



with the explicit rotation 





sin( 

cos I 




\ / 1 





(C.21) 



For example, a 1 



/I 

COS0 

— sin I 
\ 

/ 1 

cos 9 cos (f> 
— cos 6 sin i 
\ - sin 6* 

4 process can be expressed as 3 isotropic decays (see Figure KTTl) : 

^2 + ^345 ) ^345 ^3 + ^45 ; ^45 — > ^4 + ^5 ! (C.22) 



\ 

cos^ sin^ 
10 
1 / \ -sine cos^ / 

\ 

sin (f) sin 6 cos (f) 

cos — sin 6 sin 

cos^ 



where ^345 = k^ + k^ + k^, /C45 = ^4 + ^5 . The phase-space element can be written as 
d^i^4 = / d^ki, kl fc345)c/$(/c345, ^3. ^45)^^(^45, k^dkl^dkl^^ 



H.)— ri^UJ- xy\kl kl kl,,) Xy\kl,, kl kl,) X^^kl,, kl kl 



{1,2 \min 

^^45*^^345 



8kl 



( \min 

d(j)i / d cos 61 / ( 
J-i Jo 



Ofi.345 

1 ^27r 
d COS 02 



8kl 



(C.23) 



d cos 6*3 



with 



(k 



345)"™ = (^3 + Mr 



(1/2 \max 
{'^45) 



345^ 



[nil — m2 — m^y 
= {nil - ni2Y ■ 




kl /C345 /C45 /C5 

Figure C.l: The 1^4 process expressed as 3 isotropic decays 
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Similarly to the 1^4 process, the phase-space element of the 1^3 process can be 
written as, 

J (i$i^3 = J d^ki, kj, kl^)d^kl^, kl kl)dkl^ 

_ f^^i^r- xnkl kl kl,) xy\kl, kl kl) 

J{k^J"^i" 1 34 

/'27r /"l /'27r /•! 

/ d(f)i / d cos 9i / d02 / cos 02 ■ 
Jo J-1 Jo J-1 

The decay width of the 1 — > n process is defined as 

dTi^n = / . [y^lMfd^i^n, (C.25) 

where M is the matrix element of all the diagrams for this process, it is squared and 
averaged over the spin of the external particles. The production process, i.e. the 2 —>■ n 
process, is treated analogously. 
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C.3 Mult i- Channel Approach 

The contributions of the real photon bremsstrahlung can be expressed as f l4.28p . The 
singular part is separated by the phase-space-slicing method and calculated analytically. 
The finite part is calculated by the Monte Carlo method. The amplitude of the real pho- 
ton bremsstrahlung has different propagators corresponding to different diagrams. These 
propagators behave differently in different phase-space regions. Therefore, in order to ob- 
tain a stable numerical result and to reduce the Monte Carlo integration error, we use a 
multi-channel Monte Carlo method [STj 158]. 

In the case of the n-body decay, the decay width is expressed in (lC.25p . For each type 
of propagator we choose a suitable variable set Xk, the decay width can be expressed as 



dTi^n = J f{xk)p{xk)dxk, (C.26) 

where 

f(xk) = , }, , y"\M\\xk) (C.27) 

and p(xfc) is the phase-space density. Accordingly a mapping Xk = hk{r) with the random 
number < < 1 is chosen, 

dV.^n = ['^j}^dr. (C.28) 
Jo 9{hk{r)) 



The resulting density 



describes the particular behavior of this propagator. All densities g{hk{r)) are combined 
into one total density (^tot which is expected to smooth the integrand over the whole inte- 
gration region. The phase space integral flC.26p can be written as 



dTi^n 



AI „ 

/ dxkp{xk)g{xk) 

k=i 



gtot{xk 

M „i 



^Jo 9tot{hk{r)) 



The total density is defined as 



M 



gtotixk) = ^g{xk), (C.31) 



k=l 



where M is the number of the mappings (channels). 
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An example 

For the decay X2 Xi^^^^l^ we have 14 channels. Each channel smoothes a particular 
behavior of a propagator. The channels for interference contributions are not included. A 




(ki) It x\ ih) 

Figure C.2: An example of diagram of the decay X2 ~^ Xi^~^~^1 



diagram is shown as an example in Figure [(T2l The phase space integral can be decomposed 
as 



(^'35 )min 



where k 



35 



d^k,,, kl kl) j d^k2,, kl kl) , (C.32) 
^3 + /cs, ^24 = ^2 + ^4- The upper and lower limits on the variables /C35 and 



(^35)min 



(^35) cut 5 



(^35) max 
2 \max 



{rn^o — m^o — miY 



(^24)""" = {m^o + mif , (A;24j""""" = [m^o - m^^) 

where 777,35 = V^5- The infrared and collinear singularities are excluded by the cut {kl^)cnt- 
For the SPSla parameter set where li can be on shell, the variables /C35 and /C24 cire mapped 
via 

2 



k'^ 



h(ri, mf, U, {k'^^)ramj (^35)max) 5 



k'^ 



hi r 2,-1X17 

^ 1 



im 



2 \min f 1,2 \max\ 
? \'^24) I 1 



(C.33) 



where the function h has been defined in flC.6l) and flC.9p . The total density for this set of 
mappings is 

„2 2 ■ -p o /l2 \mm / t 2 \max\ „ / t 2 2 / ; 2 \ / t 2 



9 — 9x{k24, "rrtj^ im^^F^-^, 2, (/C24)™™) {^24)^'^^) 9x{k^K,, rrii , z/, {k^c^)nim, (^35)max) 



5'd(^|o, ^35, ^24)5'd(^24' ^2? ^Ds'o'l^SS' ^3' 



5^ ' 



(C.34) 



where the densities ^^a; and (7^ are defined in flC.6p . fIC.lOp and (lC.16p . respectively. The 
density in flC.34p describes the propagator of the diagram in Figure IC.2I Other diagrams 
for the decay X2 ~^ Xi^~^~^'l have been shown in Figure IT5l Their variables and mappings 
are chosen similarly to this example. 
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Appendix D 
Feynman Diagrams 



In this appendix wc present the generic Feynman diagrams for the virtual corrections of 
the decays X2 ~^ Xi^^^^i^ = s?/^?''")- These Feynman diagrams are classified into vertex 
diagrams, self-energy diagrams and box diagrams. The notations are as follows. 

Leptons = {i/e, e, i/^u, /x, i/^, r} 

Sleptons = {i>e, e, z>^, /i, z/t-, f} 

Gauge Bosons V = {7, Z, W^} 

Neutral Higgs Bosons and Goldstone Boson 0° = {h^ , H'^ , A'^ , G^} 

Higgs Bosons and Goldstone Bosons (/) = {(p'^, H^, G^} 

Neutrahnos and Charginos X — {Xi : xf} {'^ — ' " J = 1; 2) 

MSSM Fermions F = {E, u, d, c, s, t, b, x} 

MSSM Scalars 5" = (e, u, d, c, s, i,b,(f)^ 

Fadeev— Popov Ghost U — {u^,u^} . 



105 



106 



APPENDIX D. FEYNMAN DIAGRAMS 






Figure D.l: Vertex Feynman diagrams for X2 decay into Xi ^^'i two leptons via the 
the exchange of the sleptons /+, s = 1,2 labels the slepton mass eigenstates. The diagrams 
in the last line are the corresponding counterterm diagrams. 
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Figure D.2: Vertex Feynman diagrams for X2 decay into Xi ^-nd two leptons /"/"'" via the 
exchange of the sleptons Ij , s = 1,2 labels the slepton mass eigenstates. The diagrams in 
the last line are the corresponding counterterm diagrams. 
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Figure D.3: Vertex Feynman diagrams for X2 decay into x? and two leptons via the 
the exchange of the gauge bosons V. The diagrams in the last hne are the corresponding 
counterterm diagrams. 
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Figure D.4: Vertex Feynman diagrams for X2 decay into Xi ^^'^ t~t^ via the exchange 
of the neutral Higgs bosons and Goldstone boson cjp. The diagrams in the last line are 
the corresponding counterterm diagrams. The Higgs intermediate states are neglectd for 
I — e, II. 
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Figure D.5: Self-enery diagrams for the decays X2 ~^ the case of the slepton 

mixing, s{t) = 1,2 labels the slepton mass eigenstates. The diagrams in the last line are 
the corresponding counterterm diagrams. 




Figure D.6: Self-enery diagrams ior Z — V and (j)^ — V mixing. The diagrams in the last 
hne are the corresponding counterterm diagrams. 
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Figure D.7: Self-enery diagrams for 0° and Z— 0° mixing, where the Higgs intermediate 
states are neglectd for / = e, /i. The diagrams in the last hne are the corresponding 
counterterm diagrams. 
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